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INTRODUCTION 

Correlation plots provide a useful framework for analyzing several important aspects of model 
performance. The characteristics of the correlation plot between two species reflect the relative 
efficiency of transport and mixing versus photochemistry and also the relative distribution of the 
net source/sink terms for the species. The advantage of this comparison is the ability to exploit the 
high resolution sampling of normal stratospheric species fluctuations in the aircraft and balloon 
measurements. By expressing the various species as a function of a nearly conservative tracer like 
N 2 O, we can map measurements over a limited altitude and latitude range to a much broader range 
of space in the two-dimensional zonal mean representation, thereby greatly increasing the value of 
the comparison between the models and measurements. Furthermore, the balloon and aircraft data 
sets include a fairly extensive number of long-lived species with varying lifetimes and 
loss/production mechanisms, thus providing a variety of points of contact to model processes. 

In addition to the comparison of species correlations, this section also directly compares 
latitude profiles of O 3 , NO y , and simultaneously measured NOy /03 from the ER-2 aircraft with the 
corresponding quantities from the models near aircraft altitudes (about 20 km). These 
measurements represent the most extensive in situ data available for these quantities, covering 
nearly all latitudes of both hemispheres. 

The primary goal of this experiment is to provide a framework and the data with which to 
analyze model performance. The theoretical basis for interpreting correlation characteristics in 
terms of specific model processes is briefly reviewed in the next section. Theory, and includes 
guidance on the limits of this approach. The data used for this comparison are discussed under 
Measurements. In some cases, uncertainties in the measurements severely limit the value of the 
comparison to the models. In Comparison and Model Correlation the general features of the 
comparison of correlations from measurements and models are discussed and some specific 
examples are presented. Examples of the comparison of latitude profiles of NO y , O 3 , and N 0 y /03 
are presented next followed by a brief summary. 


THEORY 

A theory for the existence and significance of species correlations has been proposed by Plumb 
and Ko (1992). Stratospheric transport may be viewed as being represented by two components: 
mixing along isentropic surfaces associated primarily with breaking planetary waves and the 
planetary scale "Brewer-Dobson" circulation. The approach follows the "equilibrium slopes" 
hypothesis of Mahlman et al. (1986) and Holton (1986). A species whose stratospheric lifetime is 
long compared with the time scale of the isentropic mixing will become well mixed, not along the 
isentropic surfaces but along "rapid exchange" surfaces that dip poleward and downward relative 
to the isentropic surfaces because of the effects of the mean circulation. The location of such a 
surface can be defined by the (single) value of the mixing ratio of any long-lived reference species 
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such as N 2 O; for this value of N 2 O, the mixing ratio of any other long-lived species will also have 
a single value. Therefore, a correlation plot of these two species will be compact, i.e., the plot will 
not be scattered (though it may not be a straight line). Using a two-dimensional transport 
formulation Plumb and Ko (1992) argue that the slope of this compact curve is equal to the ratio of 
the globally integrated, vertical fluxes of the two species through the rapid exchange surfaces. In 
regions where the lifetimes of the two species are very long - compared with the time scale of 
(relatively slow) quasi-vertical transport - both fluxes will be constant and therefore the slope of 
the correlation will be linear. 

For a species to be long lived requires a lifetime of maybe a year in winter middle latitudes 
where mixing times are short, but much longer in summer and in the tropics. In the latter case, 
there may be few if any species whose chemical lifetime in the tropics is long enough to satisfy this 
criterion. 

These concepts suggest the following interpretations of the correlation diagrams. The 
compactness (or noncompactness) of the correlation indicates that the transport rates within the 
rapid exchange surfaces are short compared with the local chemical lifetime. The curvature of the 
relationship is indicative of nonconstant net flux and thus of differing quasi- vertical distributions of 
integrated sources and sinks for the two species. In correlations that are linear in the lower 
stratosphere, the value of the slope of the relationship can be used to determine the ratio of upward 
fluxes into the stratosphere and thence the ratio of species lifetimes in the stratosphere. 


MEASUREMENTS 

Table H-l lists the long-lived species considered here and their approximate stratospheric 
lifetimes, as provided by the World Meteorological Organization (WMO) (1992). The CFC 
species and N 2 O have no known losses in the troposphere, so their stratospheric lifetimes are equal 
to their total atmospheric lifetimes. These lifetimes are derived empirically and from model 
calculations. The stratospheric lifetimes for CH 4 and the hydrochlorofluorocarbons are estimated 
from model calculations (e.g., Prather and Spivakovsky, 1990; Fung et al., 1991). 

Several issues must be considered when interpreting the comparison between atmospheric 
measurements and two-dimensional steady-state, zonal mean model representations. Clearly, 
measurement uncertainties limit the degree to which agreement can be expected. The major sources 
of uncertainty in measurements of long-lived species are variations in calibration standards, 
instrument precision, and sampling efficiency. Variations in calibration standards include 
uncertainty in the absolute accuracy of the standards and possible drift of the standards with time. 
These uncertainties are manifest in systematic calibration differences between different 
measurement groups (WMO, 1992). The model uncertainties that limit the value of comparisons to 
data are generally questions of the representativeness of the model formulation. For example, 
comparison with steady-state stratospheric models is limited by the effects of annual increases in 
the tropospheric abundances of long-lived species. Systematic meridional variations, model 
boundary conditions, and parameterization of processes such as heterogeneous reactions should 
also be considered. The magnitudes of some of these uncertainties are listed in Table H-l. Their 
effects will be discussed below for each of the individual correlations. 

The data for this comparison come mainly from balloon-borne cryogenic samples (Schmidt et 
al., 1987, 1991) and whole-air samples from the NASA ER-2 stratospheric aircraft (Heidt et al., 
1989; Kawa et al., 1992). Note that the balloon and aircraft data sets are complementary in giving 
high-resolution vertical and horizontal sampling, respectively. Simultaneous measurements of 
N 2 O and CH 4 and measurements of several other long-lived chlorine species are also available 
from the ATMOS experiment (Gunson et al., 1990; Zander et al., 1987). Table H-2 lists the 
specific sampling details for these main data sources. None of the correlation data shown here are 
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from tropical latitudes. The raw data from the balloon samples, ATMOS, and the average NO^ 
O 3 , and NOy /03 from the ER-2 are available from the UADP (see Volume I, chap. 3). The aircraft 
data are on CDROM for the Airborne Antarctic Ozone Expedition (AAOE) and Airborne Arctic 
Stratospheric Expedition (AASE) (Winkler and Gaines, 1989; Hipskind and Gaines, 1990). Fits 
to preliminary data for CFCI 3 and C 2 F 3 CI 3 from the Airborne Chromatograph for Atmospheric 
Trace Species (ACATS) experiment on the ER-2 during the AASE-II are also included (J. W. 
Elkins and D. W. Fahey, personal communication, 1992). 


Table H-l. Measured Species, Stratospheric Lifetimes, and Uncertainties (from WMO, 1992) 


Species 

Stratospheric 

Lifetime 

(yr) 

Accuracy of 
Standards 
(%) 

Precision 

Balloon 

(%) 

1 

Aircraft 

(%) 

Annual 

Increase 

(%) 

n 2 o 

130 

5 

5 

5 

0.2 

ch 4 

150 

5 

1 

2 

1 

CF 2 C1 2 

120 

5 

5 

5 

4 

CFC 13 

55 

10 

5 

1 

4 

C 2 F 3 C 13 

110 

20-40* 

5-10 

1 

10 

CHF 2 C1 

240 

5-10 


15 

7 

CC14 

50 

20 

5-10 

15 

1.5 

CH 3 CC 13 

50 

10-25 

5-10 

15 

4 

CH 3 C 1 

50 

20-30 

5-10 

15 

0 

Cly 


20 # 



3 

NOy 


15# 



0 

03 


5# 



0 


* Data scaled to remove the systematic difference. 

# Estimated total uncertainty. 


As indicated in Table H-2 the measurement data span about 9 years, so that interannual 
increases in the abundance of many of the species are significant. To attain a more representative 
comparison to the models that were run for 1990 steady-state conditions, we have scaled the data 
by the amount of the annual increase for each species in the troposphere, i.e., (scaled data) = (raw 
data) (1+ fractional increase)O990-year of measurement). The fractional annual increases are derived 
from the Ozone Trends Panel and WMO repons (Watson et al., 1988; WMO, 1990,1992) and are 
shown in Table H-l. Thus, assuming that transport is similar from year to year, the scaled data 
simulate an ensemble of samples all taken with tracer values at 1990 levels. Although stratospheric 
transport is not exactly the same interannually (Schmidt and Khedim, 1991), this procedure greatly 
reduces the variability in correlations with species with significant growth rates. N 2 O values have 
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not been scaled because the annual increase is negligible compared with other uncertainties in this 
comparison. No attempt has been made to account for different ages of air in samples (at different 
altitudes, for example) for a given year. 


Table H- 2 . Measurement Sampling 



Balloon 

Aircraft 

ATMOS 

N samples 

204 

337 

19 

Year 

1982-1990 

1987, 1989 

1985 

Season 

fall, winter* 

winter 

spring 

Latitude 

44N, 68 N 

55-80 

48S, 28N 

Altitude (km) 

10-35 

16-20 

22-62 


*Also includes a few spring and summer samples 


The measurement correlations shown in this section use N 2 O as the independent variable for 
correlation plots because of its long lifetime, small annual increase, relatively simple 
photochemistry, and accurate measurements (Table H-l). A long lifetime for N 2 O has been 
corroborated by Ko et al. (1991) in comparison with Stratospheric and Mesospheric Sounder 
(SAMS) measurements although their lifetime estimate is somewhat lower (110 years) than that 
from WMO (1992). In many cases the data are fit to a line that serves as a reference for 
comparison to model correlations. The examples shown here are only a small sample of those 
possible from current data sets for models and measurements. Numerous other correlations are 
possible and may be used to examine processes specific to those species (e.g., Solomon et al. 
[1992]). 

N2O-CH4 

Figure H-l shows the correlation plot for N 2 O and CH 4 and a linear least-squares fit to the 
balloon and aircraft data for N 2 O > 20 ppbv. ATMOS data are also plotted but are not included in 
the fit. The correlation is tightly distributed and nearly linear with a correlation coefficient ( r value) 
of 0.983. The high degree of linear correlation between these species verifies that the lifetime of 
CH 4 is long in the lower stratosphere. Although N 2 O is primarily lost via photolysis and CH 4 
through reaction with OH, two processes that may have different spatial distributions, both 
lifetimes are sufficiently long in the lower stratosphere to satisfy the conditions for a compact, 
linear slope. Also, the measurement accuracy is very good. With N 2 O less than about 30 ppbv the 
correlation becomes distinctly nonlinear. Measurements in this regime are from altitudes above 
about 33 km where losses of one or both species become significant. Note that the slope of the 
correlation plot in the linear region yields a ratio of about 1.4 for the lifetime of CH 4 relative to 

n 2 o. 

N2O-CF2CI2 

Figure H-2 shows the correlation plot for N 2 O and CF 2 CI 2 (CFC-12). The data distribution is 
fairly compact, and a linear fit to the data for N 2 O > 20 ppbv gives a correlation coefficient of 
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0. 980. A second-order polynomial fit to the data (not shown) is within 10 pptv CF 2 CI 2 of the 
linear Fit for 50<N20<320 ppbv. Thus, within the regime of N 2 O values greater than about 50 
ppbv, the CF 2 CI 2 lifetime is long relative to transport, in agreement with the estimates of Table H- 

1. Values of N 2 O less than 50 ppbv are generally from samples taken above 30 km. Note that the 
CF 2 CI 2 fit to the tropospheric N 2 O mixing ratio (308 ppbv) is 407 pptv, which is significantly 
lower than measured tropospheric CF 2 CI 2 for 1990 (about 470 pptv) (WMO, 1991). This may 
reflect systematic uncertainty in the measurements; however, CF 2 CI 2 is not expected to be lost in 
the grab samples (Knapska et al., 1985). The compactness of the correlation also suggests that the 
mixing time scale is short compared with that for changing CF 2 CI 2 in the stratosphere due to the 
moderate rate of increase of CF 2 CI 2 in the troposphere. The CF 2 CI 2 values from the AASE used 
for this report are slightly different from those published on the mission CD-ROM as a result of a 
post-mission calibration comparison (L. E. Heidt, personal communication, 1992). 

N2O-CFCI3 

Figure H-3 shows the balloon and whole-air aircraft data for CFCI 3 (CFC-1 1) and the best 
polynomial fit to that data. Also shown is a fit to CFCI 3 data from the ACATS instrument (J. W. 
Elkins and D. W. Fahey, personal communication, 1992) with simultaneous N 2 O (M. 
Loewenstein, personal communication, 1992) from the ATLAS instrument (Loewenstein et al., 
1990) aboard the ER-2 during AASE-II in 1991-1992 (the CFCI 3 and N 2 0 from AASE-II are 
preliminary data subject to revision). The ACATS curve is fit to 1654 data points, and the data are 
scaled back 2 years to 1990. The ACATS data generally agree closely with the balloon and whole- 
air sample data, within the measurement uncertainties of about 10%. The variability of the data is 
comparable in both sets. All of the data show a fairly compact distribution with distinct curvature 
in the correlation relationship over nearly the entire range of N 2 O values. This compactness 
indicates that compared with mixing, the local CFCI 3 lifetime is long in the measurement regime 
while the curvature reflects the greater loss of CFCI 3 , and hence shorter lifetime, compared to 
N 2 O. The CFCI 3 loss is expected to occur mainly in the tropical lower and middle stratosphere. 
The ACATS fit (which most nearly represents in situ sampling) extrapolates more closely to the 
tropospheric measured value for CFCI 3 in 1990 (265 to 278 pptv [WMO, 1991]) at tropospheric 

n 2 o. 

N2O-C2F3CI3 

Figure H-4 shows the balloon and whole-air aircraft data, the polynomial fit to that data, and 
the fit to the ACATS/ATLAS data for C 2 F 3 CI 3 (CFC-1 1 3). The raw balloon data are analyzed 
relative to the OGIST standard (Schmidt et al., 1991), which is 37% lower than the NOAA-CNtt)L 
standard used by the whole air and ACATS groups (WMO, 1992). Consequently, we have scaled 
the balloon data by 1.37 for this comparison. The balloon and whole-air data show considerable 
scatter even after scaling to remove calibration standard differences and annual increases. Given 
the fairly long lifetime of C 2 F 3 CI 3 , this scatter is most likely due to measurement uncertainty in 
combination with the effect of the rapid annual increase of C 2 F 3 CI 3 rather than local 
photochemistry. The fit to the ACATS data is from 2034 points, which show substantially less 
scatter than that shown in Figure H-4. The fit lines are very similar, however, and both indicate 
definite curvature to the correlation relationship. Similar to CFCI 3 , this curvature reflects a shorter 
lifetime for C 2 F 3 CI 3 relative to N 2 O. Both curves extrapolate to a value somewhat lower than the 
tropospheric C 2 F 3 CI 3 value of 70 pptv for 1990 (WMO, 1992) at N 2 O = 308 ppbv, but the 
measurement uncertainties are fairly large. The rapid annual increase for C 2 F 3 CI 3 probably 
contributes to this discrepancy since the age of the lower stratospheric air is not precisely known. 

N 2 0 -CHF 2 C 1 

Data for CHF 2 CI (HCFC-22) are only available from the whole-air samples during the AASE 
and are shown in Figure H-5. A line is fit to the data to serve as a reference for comparison to the 
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models, but the fit does not explain much of the variance (r=0.58). The large scatter in the data is a 
result of measurement uncertainty derived from imprecision in the analysis (Table H-l) and the 
possible lack of conservation of CHF2CI in the sampling containers. Furthermore, CHF2CI is 
thought to have a significant interhemispheric gradient between the north and the south, while N2O 
is much more uniformly distributed. This difference may contribute to variability in the 
stratospheric measurements since air entering the stratosphere will have different amounts of 
CHF2CI for the same value of N2O depending on the hemisphere from which the air enters. The 
long stratospheric lifetime calculated for CHF2CI suggests that the variability of the measurements 
is not the result of local photochemistry. Based on Figure H-5, the CHF2CI measurements are 
only of limited value in comparison to the models. 

N2O-CH3CCI3 

Figure H-6 shows data for CH3CCI3 (methyl chloroform) and the best polynomial fit to that 
data for N2O > 90 ppbv. The data definitely indicate curvature in the correlation relation at nearly 
all N2O values, but the scatter about the fit is fairly large. The scatter in the data reflects the shorter 
lifetime of CH3CCI3 in combination with measurement uncertainties (Table H-l) and the fact that 
some CH3CCI3 is probably lost in sampling (Knapska et al., 1985). Tropospheric CH3CCI3 in 
1990 is about 140 pptv, which is significantly higher than the measurement fit at 308 ppbv N2O 
but well within the range of the data. 

N2O-CCI4 

Data for CCI4 (carbon tetrachloride) and the polynomial fit are shown in Figure H-7. Again the 
scatter about the curve is large, reflecting the difficulty of the measurement and the relatively short 
lifetime of CCI4 (Table H-l). The annual increase of CCI4 is small, and the fit value at N2O = 308 
ppbv is in good agreement with tropospheric CCI4 (108 pptv) for 1990. 

N2O-CH3CI 

Figure H-8 shows data for CH3CI (methyl chloride). The data are very scattered and suggest 
the possibility of a systematic difference between the aircraft and balloon data sets. No reference 
line has been fit for CH3CI. Tropospheric CH3CI is about 600 pptv. Although this is a major 
component of organic chlorine and the only known chlorocarbon with a significant natural source, 
the difficulty of the measurement limits its usefulness for comparison to models (Schmidt et al., 
1985; Knapska et al., 1985) 

N 2 0 -Cly 

Cly is defined as the sum of inorganic chlorine and is important to evaluating chlorine budget 
studies using measured data. The major species are HC1, CIONO2, CIO, HOC1, CI2O2, and 
atomic Cl. The origin of Cly in the stratosphere is decomposed organic chlorine molecules, the 
CFCs, whose source is near the Earth’s surface. Figure H-9a shows the correlation plot for the 
sum of measured organic chlorine (CCly) from the AAOE and AASE (Kawa et al., 1992). The 
polynomial (curved) fits have been chosen for estimating Cl y from N2O. Cl y for a given value of 
N2O is obtained by subtracting CCly for that N2O from the estimated total chlorine for the air 
parcel. Total chlorine is expected to be a conserved quantity in the stratosphere. The value for 
total chlorine is determined by the organic chlorine abundance of the tropospheric air that is 
transported into the stratosphere. Measurements of organic chlorine in the troposphere show an 
increase of about 3% per year (WMO, 1992). An average age of 3 years has been used to 
determine total chlorine for the time of the AAOE and AASE data. Figure H-9b shows the 
correlation plot for Cly based on the CCl y measurements of AAOE and AASE scaled up to 1990. 
The estimated uncertainty in the Cl y values is ±20% (Kawa et al., 1992). 
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N 2 0-NOy 

NOy is defined as the sum of reactive nitrogen species (Fahey et al., 1989). NOy is directly 
measured on the ER-2 aircraft and the principle NO^ species are all measured by ATMOS. The 
N0y-N20 correlation has been used extensively as a diagnostic for processes related to polar 
stratospheric clouds (e.g., Fahey et al., 1990). An example of the strong anticorrelation of NOy 
with N 2 O is shown in Figure H-10. The relationship is tight and linear for N 2 O >120 ppbv. The 
fits to the AAOE and AASE data, obtained from thousands of measurements during tens of flights, 
are NOy = -0.065(N 2 0)+20.6 ppbv in the Arctic and NOy = -0.074(N 2 0)+21.8 ppbv in the 
Antarctic, for (N 2 O) >150 ppbv. The ATMOS data also fit these correlations well and indicate 
that the linearity of the correlation extends down to N2O values of about 50 ppbv. The 
compactness and linearity of this correlation suggest that the lifetime of NOy is at least several 
years in the middle and high latitude lower stratosphere (in the absence of polar stratospheric 
clouds). The estimated uncertainty in the correlation slope, derived from the individual 
measurement uncertainties and the standard deviation of the fit is ±16% at the 1 c level (Fahey et 
al., 1990). 

NOy/0 3 

Latitude profiles of ER-2 measurements of NOy, O 3 , and simultaneously measured NOy /03 arc 
shown in Figure H-ll. These data are from the STEP, AAOE, and AASE measurement 
campaigns. Measurement uncertainties are ±15% for NOy and ±5% for ozone (Fahey et al., 
1989; Proffitt et al., 1989). The averages and standard deviations are calculated over 5 degree 
latitude intervals from 120-s average data points as presented in Murphy et al. (1992) with 
additional data from 870823, 870916, 870921, 8901 12, 890116, 890120, and 890124 (dates are 
read YYMMDD). Extratropical data are from the winter season in either hemisphere (except for 
870929, 871001, 871003, and 881015). Data are from potential temperatures between 430 and 
520 K, which corresponds to 45 to 65 mb in the tropics and 50 to 80 mb at middle and high 
latitudes. The advantage of using NOy/03 for comparison of horizontal gradients can be seen in 
the lower variability of NOy/03 relative to the mean which is lower than that for either species 
individually. With respect to the correlation figures shown above, the NOy/0 3 plots represent a cut 
across latitude of the slope of the NOy-03 correlation in the lower stratosphere. The ratio will have 
variations with latitude (approximately isentropically) where the local production/loss processes of 
these species become significant relative to mixing transport. For further discussion see Murphy et 
al. (1992). 


COMPARISON WITH MODEL CORRELATIONS 

The model scenario specified for comparison of tracer correlations is shown in Table H-3. 
Note that the tropospheric boundary conditions for the tracer mixing ratios are fixed at values 
representative of approximately 2.5 years prior to 1990. Since the measurements are scaled by 
annual increases to be representative of 1990, measured values for air with a residence time of 2.5 
years in the stratosphere are expected to give the best agreement with the steady-state model 
results. These should be values in the middle of the N2O range of the correlation plots. Species 
with small annual increases are expected to give closer comparisons than those with large 
increases. 

The model groups participating in this comparison are CAMED, DUPONT, GSFC, ITALY, 
LLNL, MRI, NCAR, and WASH; AER and NOCAR provided a partial set of results. Model 
abbreviations are defined in chapter 3. The large number of plots generated for this comparison 
cannot totally be accommodated in the final report. We have tried to select a representative subset 
of examples to demonstrate the features discussed below. 
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Table H-3. Bulk Tropospheric Mixing Ratios for 'Current' Stratosphere Runs 


yearSS 

Fll F12 FI 13 FI 14 FI 15 1211 1301 H22 CCI 4 CH 3 CCI 3 

CO 2 

n 2 o 

CH 4 

1980SS 

1990SS 

149 250 11 3 

253 434 44 7 

2 0.4 0.4 41 93 85 

5 2.0 2.6 92 103 145 

— > ppt 

334 

350 

ppm 

300 

308 

PPb 

1538 

1685 

PPb 


PPl- 


also fixed: CH3C1=600 ppt. 

CH3Br=15ppt, N2=78%, 02 = 21 % 




These steady-state scenarios use tropospheric boundary conditions from 2.5 years previous and are 
identical to the UNEP scenarios. Total Cl=3.3 ppbv 


General Features 

Most general features of the model results are consistent with expectations for long-lived 
species. Very slight seasonal differences are found between the model correlations from December 
and March, and differences between gas phase and heterogeneous chemistry results are negligible 
(including species lost in reaction with OH). Also, the hemispheric asymmetry is small for most 
comparisons; the primary latitude dependence is the difference of the tropical points from those at 
higher latitudes, consistent with expectations from theory outlined above and by Plumb and Ko 
(1992). These general features can be observed in Figure H-12. Because of the similarity of the 
March and December models and also the gas phase and heterogeneous cases of each model, only 
correlations for December, gas phase simulations will be shown in the comparisons that follow. 

When comparing the data for most cases, the differences between models are not much greater 
than the variability in the data. Although differences exist, the data are often not precise enough to 
indicate which models are more accurate. The models differ substantially in the compactness of 
their tropical points, but little tropical data are available for comparison. All the models seem to 
have weak poleward-downward circulation at high latitudes compared with the data. Steady-state 
models are not very useful for comparison to species with significant growth rates. These features 
are consistent with the results of the model-model comparison of idealized tracers in section O. 

CH 4 

All of the models do reasonably well in simulating the N2O-CH4 correlation in absolute value 
and slope in the linear region of N2O greater than about 50 ppbv (Figure H- 13 ). The agreement is 
consistent with the long stratospheric lifetime and small annual rate of increase for CH4 (Table H- 
1). The models differ chiefly in the amount of scatter in the correlation but, except for the tropical 
points, almost all the model points are within the envelope of the measured data points. The 
tropical points are expected to form a more scattered correlation as discussed above. The models 
also appear to do well in simulating the nonlinear decrease of CH4 at N2O less than 50 ppbv in 
comparison with the few measured points. Comparison of N2O and CH4 latitude-height cross 
sections and vertical profiles along with correlation plots from SAMS data are found in section E. 

CF2CI2 

With reference to the N 2 O-CF 2 CI 2 correlation, the models all differ from the measurements in 
predicting higher CF 2 CI 2 values for N 2 O greater than about 200 ppbv (Figure H-14). This 
difference follows from the difference between the extrapolation of the measured correlation and 
CF 2 CI 2 values measured in the troposphere as discussed in the measurement section above. The 
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models all predict a relatively tight, nearly linear relationship similar to the measurements, as 
expected for the long lifetime of CF2CI2. 

CFCI3 

The model correlations for CFCI3 exhibit significant curvature and scatter, in general 
agreement with the measurements, with the tropical points deviating furthest from the reference 
lines (Figure H- 15 ). The models are all similar in this respect except for ITALY, which has most 
points below the reference lines and greater scatter, indicative of slower horizontal transport 
relative to the photochemical CFCI3 sink. The relatively high amount of scatter is a consistent 
feature of the ITALY model. 

C2F3CI3 

The comparison for C2F3CI3 is dominated by the effect of the rapid annual increase in the 
troposphere for this species (e.g., Figure H- 16 ). The model values all lie below the reference 
curve for N2O greater than about 100 ppbv, because their steady state tropospheric boundary 
condition of 44 pptv is much less than that of ambient C2F3CI3 in 1990 , i.e., 70 pptv. The model 
correlation slopes are all less steep than the measurements, as expected for comparison of steady- 
state versus increasing tropospheric abundances. C2F3CI3 is clearly not a good point of 
comparison for models using steady state tropospheric input, but could be of considerable value in 
testing models with varying tropospheric input. 

CHF 2 C1 

The correlation comparison for CHF2CI is not constrained very closely because of the large 
scatter in the measurements, much of which is attributed to measurement uncertainty. Among the 
models, however, the differences for this correlation are large. Figure H -17 shows the range of 
values and scatter for the models. These differences indicate significant differences in the models’ 
transport versus loss due to OH reaction, the primary loss for CHF2CI. Such large differences do 
not appear in the correlations for CH4 or CH3CCI3, other species lost in reaction with OH. 
CHF2CI is an important CFC replacement compound, so a better measurement data base and 
careful analysis of model loss processes for this compound are very important. 

CH3CCI3 

CH3CCI3 is lost both through photolysis and reaction with OH but its lifetime is much shorter 
than CHF2CI. This leads to greater scatter in the correlation relationship for CH3CCI3 compared 
with CHF2CI in the models, but in this case the model results are all rather similar to one another 
(e.g., Figure H- 18 ), more like CFCI3. Most model values (outside the tropics) are higher than the 
reference curve, consistent with the possible sampling loss of CH3CCI3 in the measurements; 
however, comparison with Figure H-6 shows that the model points fall within the range of the 
measured data. 

CCI4 

The correlation for CCI4 shows a large amount of scatter in all the models (e.g., Figure H- 19 ), 
especially in the tropical points, which is indicative of its relatively short photolytic lifetime. The 
measurements are also highly scattered, but even more so than the models, due in part to 
measurement errors. CCI4 is increasing slowly in the atmosphere but the model slopes are all 
steeper than the measurements, with most points below the reference line. This difference is 
probably not significant in light of the measurement uncertainties. 
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ch 3 ci 

Comparison with CH3CI is difficult because of the large measurement uncertainty. The models 
all show a very tight correlation relationship (e.g.. Figure H-20) in contrast to the highly scattered 
measurements. There is a significant degree of curvature, indicating a difference in the 
distributions of loss rates for N2O photolysis versus reaction of CH3CI with OH. The degree of 
curvature differs somewhat among the models but not as much as for CHF2CI. 

Cly 

The model results show generally good agreement compared with Cl y derived from 
measurements (Figure H-21). Some models have slightly more curvature than the reference line 
and differ in the amount of scatter in the correlation, but most extratropical points fall within the 
scatter of the measurements (Figure H-9). The reason that the CAMED model at low N 2 0 does 
not approach the total chlorine amount specified by the boundary conditions (Table H-3) is not 
presendy apparent. 

NOy 

Most models simulate the N20-N0y correlation fairly well (Figure H-22); the CAMED model 
does exceptionally well. All the models show a nearly linear, fairly compact relationship for N2O 
greater than 100 ppbv, away from the tropics. However, several have higher NOy values than the 
reference, suggesting a problem in simulating the NOy sources and sinks. In the NCAR model, 
for example, the net amount of NOy entering the stratosphere from the troposphere appears to be 
too large. Another possible reason for overestimation of NOy in the models is loss of NOy in the 
stratosphere due to sedimentation of polar stratospheric cloud particles; this process is not 
represented in the models. The global magnitude of NOy loss on cloud particles is not well 
known. The models all show a change in slope at N2O less than about 50 ppbv, a feature also 
seen in the ATMOS measurements at high altitudes. At altitudes above about 30 km, the loss of 
NOy from both the photolysis of NO and the reaction of NO with N becomes significantly rapid 
(Fahey et al., 1990). 

NOy, 0 3 , and NOy/0 3 LATITUDE PROFILES 

Among the models a wide range of differences can be seen with respect to NOy, O3, and 
simultaneously measured NO y / 0 3 as functions of latitude (Figure H-23). For NOy and O3, the 
upper set of three model curves (20 km) corresponds more closely to the data in tropical latitudes 
while the lower set (18 km) more closely corresponds to the high latitude data. For NOy/ 0 3 , the 
lower set of model curves in the tropics is at 20 km. The December to March curves are most 
representative of the data in the northern hemisphere, while the September curves are more 
representative of the southern hemisphere data. 

Within the large differences between models some general features appear. In contrast to the 
correlation plots, seasonal variations and hemispheric asymmetries are observable in most models. 
These variations arise mainly from O3 changes, as expected from the greater photochemical activity 
of O3 relative to the longer-lived species. Also, differences are noticeable between the gas phase 
and heterogeneous simulations, especially in O3. For most models these differences are small, 
with slightly increased O3 at high latitudes in the heterogeneous runs (not shown) and a negligible 
change in NOy (DUPONT is an exception showing removal of NOy and a larger increase in O3 
than the others, with a corresponding decrease in NOy/ 0 3 , in the heterogeneous case). 

All of the models, except DUPONT gas phase, tend to predict NOy mixing ratios that are too 
low at high latitudes in comparison to the measurements (see also the comparison to LIMS data in 
section F on NOy absolute stratospheric abundance and distribution). Presumably, this feature is 
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due to the models’ weakness in simulating the polar vortex. Low modeled NO v at high latitudes is 
consistent with comparisons for other tracers in zonal mean cross sections which show that the 
poleward downward slope of the isolines is not steep enough, a feature that seems to indicate that 
the model meridional circulation versus lsentropic mixing is not strong enough. Comparison to O 3 
data at high latitudes is closer for most models, as it is at all latitudes. This leads to 
underestimation of NOy /03 at high latitudes for most models. 

m -HHi he i m ^ 1S ge n?u r ^ y A^i t0 T, u * ate wel1 the stee P N °y°3 gradient between the tropics and 
middle latitudes. The ITALY model is unique in simulating this gradient, consistent with weaker 

horizontal diffusion which seems to be characteristic of the ITALY model. However the 
magnitudes of the NOy/C >3 values from the ITALY model are too low at tropical and subtropical 
latitudes. Many of the models underpredict NO/O3 at middle latitudes. In the tropics all models 
either are in agreement with NOJO3 or overestimate the observed values (mostly due to 
overestimated NOy), except ITALY as mentioned and MRI which underestimates NO v almost 
everywhere. The general difficulty in simulating the NOy, O3, and NOJO3 data with latitude, and 
the wide disagreement between models, indicates that the two-dimensional models have serious 
problems with transport in the lower stratosphere. 


SUMMARY 

We have attempted to provide a method to analyze model performance in comparison with 
measured data. The theoretical framework and observational data have been put forth along with 
some example comparisons. We expect that the main value of this section will come in detailed 
analysis by the individual modeling groups studying their model characteristics compared with the 
measurements. In this way improved model representations of transport and chemical 
loss/production will enhance confidence in the model assessments of the atmospheric effects of 
perturbations such as stratospheric aircraft. 


An equally important aspect of this comparison has been identification of the need for 
improvement of the measurement data base. In some cases the measurement uncertainties are too 
large to rigorously constrain the model results. This situation should be improved by 
mtercompanson of calibration standards among the measurement groups and new, more precise 
mstniment techniques. Additional data, especially in the tropical regions where almost no data are 
available, will also increase the value of these comparisons. 
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Figure Captions 

Figure H-l. Correlation plot for N2O and CH4 and a linear least-squares fit to the data for N2O 
> 20 ppbv. The equation of the reference line is CH4 = 3.82(10-3)(N2O) + 0.5 ppmv where N2O 
is in units of ppbv. The standard deviations of the slope and intercept are 3.12(10-5) and 5.9(103), 
respectively. 

Figure H-2. Correlation plot for N2O and CF2CI2. The equation of the reference line is CF2CI2 
= 1.5(N20) - 55.5 pptv where N2O is in units of ppbv. The standard deviations of the slope and 
intercept are 0.013 and 2.4, respectively. 

Figure H-3. Correlation plot for N2O and CFCI3. The equation of the reference curve for the 
balloon and whole-air data is CFCI3 = 3.36(1 0-3)(N2O) 2 - 0.218(N2O) - 10.4 pptv where N2O is 
in units of ppbv. The standard deviations of the coefficients are 1.6(10-4), 0.06, and 5.8, 
respectively. 

Figure H-4. Correlation plot for N2O and C2F3CI3. The balloon C2F3CI3 data are multiplied by 
1.37 to account for the different calibration standards used by the different measurement groups. 
The equation of the reference curve for the balloon and whole-air data is C2F3Q3 = 3.93(104 XN2O) 2 
+ 0.10(N2O) - 5.25 pptv where N2O is in units of ppbv. The standard deviations of the 
coefficients are 4.9(10-5), 0.02, and 1.7, respectively. 

Figure H-5. Correlation plot for N2O and CHF2CI. The equation of the reference line is 
CHF2CI = 0.166(N2O) + 42.3 pptv where N2O is in units of ppbv. The standard deviations of the 
slope and intercept are 0.02 and 2.94, respectively. 

Figure H-6. Correlation plot for N2O and CH3CCI3. The equation of the reference curve is 
CH3CCI3 = 2.87(10-3)(N 2 O) 2 - 0.63(N 2 O) + 35.3 pptv for N 2 0 > 90 ppbv. The standard 
deviations of the coefficients are 1.4(10-4), 0.06, and 5.4, respectively. 

Figure H-7. Correlation plot for N2O and CCI4. The equation of the reference curve is CCI4 = 
1.02(10'3)(N20) 2 + 0.054(N20) - 15.0 pptv where N2O is in units of ppbv. The standard 
deviations of the coefficients are 1.4(10-4), 0.05, and 5.0, respectively. 

Figure H-8. Correlation plot for N2O and CH3CI. 

Figure H-9. (a) The organic chlorine (CCly) sum versus N2O from whole-air sample 

measurements in the Arctic and the Antarctic from Kawa et al. (1992). The solid and dashed 
curves are fit to the data by regression for a second degree polynomial and a line, respectively. 
The equations of the polynomial fits used to predict CCly from continuous N2O measurements are 
(CCly) = 2.79(10' 5 ) (N2O) 2 + 1.73(10-3) (N2O) + 0.137 in the Arctic and (CCly) = 3.19(10- 5 ) 
(N2O) 2 + 0.76(10-3) (N 2 0) + 0.013 in the Antarctic, where (CCly) and (N2O) are in ppbv. 
(b) Correlation plot for N2O and Cly derived from the measurements of CCly and estimates for 
total chlorine appropriate to the measurement regime. The equations of the reference curves are Cly 
= 3.3-(2.98(10-5)(N2O) 2 + 1.65(10-3)(N2O) + 0.13) pptv for the AASE and 3.3-(3.57(10- 
5)(N 2 0) 2 + 0.85(10-3)(N 2 0) + 0.02) pptv for the AAOE. 

Figure H-10. Example correlation plot for N2O and NOy from the flight of 12 January 1989 in 
the AASE. Data are 10-s averages from a latitude range of 59 to 79° and a potential temperature of 
460 ±10 K. From Fahey et al. (1990). 

Figure H-ll. Latitude profiles of (a) O3, (b) NOy , and simultaneously measured (c) NOy/03 
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from the NASA ER-2 aircraft. Points are 120-s averages of 1-s data. Averages (triangles) and 
standard deviations (bars) are calculated over 5 ° latitude intervals. 

Figure H-12. N 2 O versus CH 4 for the DUPONT model. Model results are interpolated onto the 
standard UADP grid and species values are plotted at every second latitude grid point for each 
altitude level. The different symbols represent different latitude ranges. Triangles are for 90 to 
30S, squares are for 20S to 20N, and circles are for 30 to 90N. Correlations at approximately the 
solstice and equinox are shown. The plots labeled het are for the model including heterogeneous 
reactions of N 2 O 5 and CIONO 2 with H 2 O on background sulfate aerosol particles. The reference 
line is a fit to the measured data as discussed in the measurements section above. 

Figure H-13. N 2 O versus CH 4 for the December gas phase simulation of each model as in 
Figure H-12. N 2 O values greater than 308 ppbv in the GSFC model are a numerical artifact 
occurring near the north pole. 

Figure H-14. N 2 O versus CF 2 CI 2 (CFC- 12 ) for each model as in Figure H-13. 

Figure H-15. N 2 O versus CFCI 3 (CFC- 1 1) for each model as in Figure H-13. 

Figure H-16. An example of N 2 0 versus C 2 F 3 CI 3 (CFC-1 13) from the CAMED model as in 
Figure H- 1 3. 

Figure H-17. N 2 0 versus CHF 2 CI (HCFC-22) for each model as in Figure H-13. 

Figure H-18. An example of N 2 O versus CH 3 CCI 3 from the LLNL model as in Figure H-13. 
Figure H-19. An example of N 2 O versus CCI 4 from the NCAR model as in Figure H-13. 

Figure H-20. An example of N 2 O versus CH 3 CI from the NOCAR model for January. 

Figure H-21. N 2 O versus Cl y for each model as in Figure H-13. 

Figure H-22. N 2 O versus NO y for each model as in Figure H-13. 

Figure H-23. O 3 , NO y , and simultaneously sampled N 0 y /03 as functions of latitude for the 
contributing models. Circles with error bars are measured averages as discussed in the 
measurements section. Lines are model results from March (dotted), September (dashed), and 
December (solid). Upper set of 3 lines for O 3 and NCk, is at z* = 20 km and lower set is at z* = 
18 km. The measurements correspond most closely to the lower dashed curve south of about 30S, 
to the upper solid curve between about 30S and 30N, and between the lower solid and dotted 
curves north of about 30N. For NO y /C >3 the lower set of lines in the tropics is for 20 km and the 
curves should bracket the observations. 
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INTRODUCTION AND BACKGROUND 

The goals of the carbon- 14 and strontium-90 radionuclide experiments are to conduct 
fundamental tests of dynamical transport in the models that are independent of chemistry. 
Previous studies have used carbon- 14 (Johnston et al., 1976; Kinnison, 1989; Shia et al., 1989; 
Jackman et al., 1991) and strontium-90 (Kinnison, 1989) as tracers in one- and two-dimensional 
models. These tracers were produced by the aboveground nuclear tests of the late 1950s and 
early 1960s. Carbon- 14 is produced by fission and fusion nuclear processes, while strontium-90 
is produced only by fission processes. The atmospheric nuclear test moratorium between the 
United States and the U.S.S.R. ended aboveground nuclear tests in December 1962. After 
January 1, 1963, the distribution of these tracers became dependent only on the transport 
processes in the stratosphere. After June 1967, the French and Chinese performed a small 
number of atmospheric nuclear detonations in which the bomb-rise height was primarily below 
20 km altitude. The injection of these radioactive species into the lower stratosphere is 
somewhat similar to the proposed high-speed civil transport (HSCT) emissions, and the time 
scales for removal of these radionuclides is a suitable test of the models. 

This section is divided into two parts: Modeling Carbon- 14, which is in the form of CO 2 and 
acts as a passive gaseous tracer; and Modeling Strontium-90, which sticks to aerosols and 
therefore must include a settling velocity to calculate its spatial and temporal variation. 
Throughout this section the comparison of model-derived tracers is compared to observed 
distributions. All altitudes reported in this chapter are pressure altitudes (z*-coordinates). In 
addition, model intercomparisons are examined. 


MODELING CARBON- 14 

Extensive measurements of excess carbon- 14 from nuclear weapons testing were conducted 
aboard U.S. aircraft and balloonsondes in the troposphere and stratosphere (see Hagemann et al., 
1965, 1966; Telegadas, 1971; Telegadas et al., 1972) as a function of time, latitude, and altitude. 
Typically, there were four latitudes (70N, 3 IN, 9N, and 42S) where carbon- 14 data were 
measured. Measurements were observed at altitudes greater than 20 km only at latitude 3 1 N. 

Telegadas (1971) reported the measured, atmospheric, excess carbon-14 data from March 
1955 through July 1969. Using these data, his experience with other radioactive isotopes (Sr, Cs, 
Zr, etc.), and his aircraft-based observations of early stage bomb clouds, Telegadas constructed 
contour plots of carbon- 14 mixing ratios (that is, units proportional to mixing ratios) for each 
quarter year over the 15 years. Each plot represents an average of all aircraft and balloon data 
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for a 3-month period (e.g., October is an average over September, October, and November). In 
Figure I- la, the altitude-latitude distribution of carbon- 14 for October 1963 is shown, including 
numerical values of measurements at four latitudes, crosses where measurements were made at 
other latitudes, and contour lines drawn by Telegadas. Johnston (1989) concluded that October 
1963 provided good initial conditions for testing the dynamical representation of two- 
dimensional models and used Telegadas' contour plot for this purpose. October 1963 was 9 
months after the conclusion of the nuclear bomb test series and was the first time at 3 IN that the 
local carbon- 14 mixing ratios at all measured altitudes passed through their maximum value. 
Since the nuclear tests were in the northern hemisphere and in the tropics, there were large 
spatial gradients of carbon- 14 at that time. Values outside the range of Telegadas’ contour lines 
for October 1963 were extrapolated by Johnston (1989), largely on the basis of observed trends 
during the period 1958 to 1961. Using the observed carbon-14 distribution for January 1960 as 
initial conditions and two models for nuclear bomb-cloud rise, Kinnison (1989) calculated the 
carbon-14 distribution as a function of time during and after the 1961-1962 nuclear bomb-test 
period having each bomb contribute at its appropriate latitude and time. The model results using 
the bomb-rise model of Seitz et al. (1968) agreed very well with Johnston's initial conditions for 
October 1963 at 70N, agreed fairly well at 31N, but gave substantially larger carbon-14 mixing 
ratios in the tropics and southern hemisphere, especially at altitudes above 25 km. On this basis 
Johnston (1989) prepared a second set of extrapolated values (Figure I-lb), having larger carbon- 
14 mixing ratios than the first set but still having values substantially lower than Kinmson's 
(Figure I-lc) in the region of extrapolation. For this study, Johnston's higher values (Figure I-lb) 
were used at altitudes above Telegadas’ contour lines. Inspection of Figure I- la shows that the 
data for the initial conditions (October 1963) are extremely sparse, but they are reinforced to 
some extent by Telegadas' experience and by Kinnison's modeling (1989). 

The time-dependent lower boundary values specified by the models for the northern and 
southern hemispheres were obtained from Johnston (1989) and are given by the following 
functions: 

14c (N-Hem.) = 73.0 - 0.27823 t - 3.45648x10-3 t2 + 4.21 159x10-5 t3 

14c (S-Hem.) = 44.5 + 1.02535 t - 2.13565x10-2 12 + 8.61853x10-5 t3 

where t is the time in months after 15 October 1963. These functions are based on observed data 
during this time period. For all times after 15 June 1968, the lower boundary values for the 
southern hemisphere were set equal to the calculated values from the northern hemisphere. The 
upper boundary conditions were specified as zero flux. The observed distributions are reported 
in units of 105 atoms of excess carbon-14 per gram of dry air, which is proportional to the mixing 
ratio, and these units are referred to as "mixing ratio units" or at times simply as mixing ratios. 

The models were used to derive the distribution of carbon-14 between October 1963 and 
January 1971. The three-dimensional GISS model was integrated only from October 1963 to 
October 1968. Comparisons are made between observed profiles at the four latitudes during the 
period from 15 October 1963 to 15 July 1966. Recall that each observed profile is an average of 
numerous aircraft and balloonsondes over a 3-month period. The variability of carbon- 14 in a 
selected number of observed profiles at 3 IN is represented in Figure 1-2. Four different times are 
included in this plot; the individual measurements along with the profile derived by Johnston 
(1989) from Telegadas' contours are shown. There is low variability in the troposphere but 
significant variability in the stratosphere. 

Observed carbon- 14 profiles are available (Johnston, 1989) for the four latitudes every 3 
months through July 1966. There are additional profiles between July 1966 and July 1969, 
which have no data above 20 km altitude and are not discussed in this section. Figures 1-3, 1-4, 1- 
5, and 1-6 give observed and calculated carbon- 14 mixing ratio vertical profiles, respectively, at 
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31b J’ ? N ’ and J 42S - In Fi g ure I *7, altitude-latitude contour plots of carbon- 14 are shown for 
each model one and two years after the initial condition (October 1964 and October 19651 A 

S—I campaign with balloon flights above 20 km at five latitudes was carried out in 
1970-71 (Figures 1-8 and 1-9) (Telegadas et al„ 1972). 

Inter-model Comparisons 

Global Stratospheric Residence Times 

„ 8 i 1 °K a i 1 invento ^ excess carbon- 14 molecules was calculated for each 3-month period 

In 8 ™ S ,T,! Phe r C Kr denCC u mes were found by a 'east-squares fit of Infinventory) vs 
time. In Table 1-1, the global stratospheric residence times (between 16 and 60 km) are shown 

“IT inle T vals - In Table 1-2, for each interval, the models are grouped according to 
residence time: short, medium, and long. Between October 1963 and 

Scth^ anH^hP^nq I i a) ’ l 5 ere are eight models in the shon stratospheric residence time 
section and the range of residence times among all the models is 1.3 to 3.2 years. When the 

ioA/^ S »K henC ^denee time for the same models are compared between January 1965 and July 

Hmt 6, ;liH e ?h m0de S ~r GIS n’ < ? SFC ’ and N CAR— have a relatively short stratospheric residence 

October i 963 Ind^.l^iQ^ 6 h ^ ' Vi? t0 4,9 . years ' In Table F2c > the period between 
VPQfC / i 1 y 1966 1 ^ shown » and range in stratospheric residence times is from 2 3 

Iannarl G /Q7 S i !°hM, (ITALY)> Table F2d s P ans the full time, October 1963 through 
1Q^ 9 IL’ T nd the res * d ® nc ^ times range from 3.8 to 4.9 years. For the time interval January 
965 though 3a, JP a r ^ ' 1 97 1 , the range of calculated residence times is 4.3 to 5.4 years. In 

D , U , PONT ’ and WASH have the longest stratospheric lifetimes. LLNLND and 
WASH have short lifetimes over the first year but show considerably longer carbon- 14 removal 
imes at longer times after the October 1963 initial conditions. As can be seen from the figures 
the models remove carbon-14 relatively rapidly from the lower stratosphere, but the middle to 
upper stratospheric removal rates are much slower. 

Altitude-Latitude Contour Maps 

Figure I- lb gives the altitude-latitude contour plot of carbon- 14 mixing ratio units used by all 

^rnhtpH 1 h nitl f C °^r/ f0 A° C k 0be , r nL 9 , 63- Figures F7 g ive altitude-latitude contour plots 
calculated by 12 models for October 1964 and for October 1965. A rough measure of the 

agreement and disagreement among the models for October 1964 was carried out bv 

^ P ! n ? ipOSing S meS f ° r C ° nt0Ur 200) Contour 300 ’ contour 40 °, and cont °ur 500, each on a 
SSTtSni S e . s t u P e H nm P os f d lines for contour 200 spread from pole to pole with maximum 

ofahnnt 4^ Thf^m es ’ and at each latltude there 1S a vertical spread, among the 12 models, 
j' Tbe , contours appear at both low altitudes and at high altitudes. The low 
altitude band was similar in shape to the band of 200 contours, and its vertical spread was about 
5 I?!'- I, ^PP er altl J tude 300 contours appear as a random tangle from 30 to 60 km in middle 
and high latitudes and from 40 to 60 km in tropical latitudes. The superimposed 400 contours 
stretch from the North Pole to southern midlatitudes, show a lower altitude band of about 7 km 
vertical spread and an upper altitude band of about 15 km vertical spread. One model has a 400 
contour line as high as 50 km and another has no 400 contour line above 30 km The 500 
contours more or less uniformly fill an envelope with a 15 to 20 km vertical spread and 

fr0I H-f? e N ° rth Po,e to about 15S - This rou i h preview of the model calculations 
indicates large differences among their results. 
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Comparisons Between Calculated and Observed Excess Carbon- 14 
Lower Stratosphere 

The boundary values at the Earth's surface are the time-dependent observed carbon- 14, which 
implicitly gives a flux into the surface, and there is no significance to the agreement between 
models and data in the troposphere. The altitude at which carbon- 14 begins sharp increase is 
a measure of the tropopause height, both for the data and for the models. At 70N, 9N, and 2b, 
the measurements were made only from aircraft, and the maximum altitude was 20 km. lhe 
observed data at 70N, Figure 1-3, appear in some cases to reach the maximum carbon- 14 mixing 
ratio and give a good comparison of models and data in the lowest stratosphere, but the 
tropopause is so high at 9N and 42S that the comparison between models and data is of limited 
value at these latitudes (Figures 1-5 and 1-6). From about April 1964 to April 1965, the altitude 
of the calculated carbon- 14 tropopause is lower ( - ), about the same ( 0 ), or higher ( + ) than that 
observed, as follows: 

Eat. Model LaL MMfii 


7 ON NCAR 

CAMED - 
DUPONT 0 
ITALY 0 
AER 0 
LLNL 0 
WASH + 
GSFC + 
CALJPL + 
GISS + 
LLNLND + 
MPI + 


3 IN ITALY - 
CAMED - 
DUPONT 0 
LLNLND 0 
AER 0 
LLNL 0 
WASH + 
GSFC + 
CALJPL + 
GISS + 
NCAR + 
MPI + 


A quantitative comparison of the data and the models is given by Table 1-3, where the ratio of 
calculated carbon- 14 to observed carbon- 14 at 3 IN latitude and 20 km altitude (the altitude of 
the observed maximum carbon- 14 mixing ratio) is given from January 1964 through October 
1965. On October 1963, all models agreed with each other and with the data, and at January 
1964 three months later, most models have a calculated/observed ratio slightly below one, with 
the range of values being 0.57 to 1.00. During the next 2 years, this ratio decreases to about 0.35 
for MPI; 0.5 for GISS, NCAR, and GSFC; 0.6 for CALJPL and WASH; and 0.7 for AER, 
LLNL and DUPONT. For three models the ratio, calculated/observed, remained fairly constant 
for 2 years: LLNLND averaging 0.90, CAMED averaging 0.97, and ITALY averaging 1.25. For 
all the models, the average ratio varies from 0.36 to 1.23. A similar comparison is given for 70N 
latitude and 16 km altitude in Table 1-4, but these data show greater irregularities than those at 20 
km. GISS, CALJPL, and GSFC decrease to about 0.6. Within the large scatter, AER, LLNLND, 
and WASH show a tendency for the ratio to decrease with time; the ratio averages about 0.9 for 
these three models. The ratio, calculated/observed, averages about 1.1 for NCAR, CAMED, and 
LLNL and it is 1 2 for ITALY and 1.4 for DUPONT. For all the models, the average ratio 
varies from 0.56 to 1.4. At (16 km, 70N) and at (20 km, 31N), the 2-year average model results 
are higher than the data for some models and lower for others. The model predictions, involving 
only atmospheric motions and no chemistry, show a factor of two spread among the 2-year 
average, lower stratosphere, sweep-out times (Tables 1-3 and 1-4). 
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Middle Stratosphere 


A ‘ 3 ' N latitude, the observed carbon-14 data show a bell-shaped vertical profile skewed 
IQdd? 'os r niCa - whose alidade of maximum mixing ratio increased from 2 lkm in Januarv 
dfffrmm h n™ Janu ^ 1966 < Tablc •-«. The model-calculated vertical profiled show qute 
^havtor; in all cases the altitude of maximum carbon-14 mixing ratio increased from 

(^t^ Ct i'l 9 ^ l ” X ^^ l ^ S ^ ^ n ^ s3 d^^i'li^^Piefdes 1 in^Figurtf M^for example ^Apri/ 1 964 S 

a^ve 

Phetr e 1aS ltS maK,mum at 23 km. By July 1965, all models but LLNLND and ITALY have 
m^ira'X^C"" 13 ^ 4 ° km ' a " d * A P"' 19 “ — '»» -odelsIlXve !S 

Another representation of this effect is shown in Table T fS Thp ttat v i * i 

rCAuVLm^?"' 14 W,h re r S r- *° 

5f CAUBLmode is an average of the other models in this respect (Table 1-5) The slnnenf 

obse^edSma be,W “" 25 a " d 32 ^ are preslnmd fP Ae 

IQttTh?,’ h I J, Y j T1 e del : and the CALJPL model in Table 1-6 from January 1964 to Julv 

in^Januarv S ^ ^ obs ?™ ed data 3,1 have a decrease of 32 mixing ratm units per km 

in January 1964, and a year later the carbon-14 data have the same vertical slnnc nnrifrr 

r^rrjL 196 ! 4 '?“• lhe “PPe ra ™ .of 

of if units to -18 to -1 1 to -2, which is almost straight up. The ITALY model 

^^H S i r ^ SU * S tbat are intermediate between the observed data and the results of the CAI IPI 

U S'n? C,0ber 1964 a " d rotatet^upwwdly^pa^ing^the 

observado'ns ^^all^of^h^modelsP Tn^lhe^'auPhcms^f^dbs'^ction 'gWe^^discussi^^of'hs 
posstble cause. These three suggestions are given as the next thn£ p^graph s d,SCUSS10n of ,ts 

t , Sfl ieWei u A ; With , the initial conditions a11 models used, there is no way the values above 
l^erta^ues b abofe a the°^ neaV^f, ? b , md ~ no matter what K " Profile one uses. One could get 
conditions are such that the equaTorial^Tluesare Iwenha^ Thes^ 

dynam™ asdfagnLS "by ’SSaT™? may E' 

be important to an understanding of Figure I- 10." y 11 


ht,cJ| he SOUrCe of disagreement between the observed carbon-14 data and all the models mav be 
based, in part on too fast transfer to the troposphere by the models as can be seen fmm 

equjva"em°to CALJPL)""* miXi " g “ 3IN inS ‘ ead of ‘ he ra ' ios <>f ^eM ^ 


Date 4/64 

Obs (20 km) 580 
JPL (20 km) 560 
Obs (32 km) 325 
JPL (32 km) 450 
JPL (40 km) 350 


7/64 

10/64 

1/65 

520 

500 

520 

380 

330 

300 

340 

300 

400 

520 

500 

460 

410 

460 

440 


4/65 

7/65 

10/65 

400 

340 

340 

280 

200 

190 

340 

320 

340 

420 

390 

360 

420 

410 

390 


1/66 

4/66 

7/66 

300 

300 

280 

200 

185 

150 

330 

300 

280 

325 

305 

290 

360 

340 

305 
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A n i, m oititnrlp the observed mixing ratios are essentially constant at 350±50 from April 
1964 to April by W 

larger that, the observed values 1^ k™ km altitude, the calculated 
and decrease during the next^ years to , y October 1964, and then decrease 

results at various times is as follows. 


9N 

42S 


Talc ;> Q1 e 

1/64 

1/64 - 1/65 


Calc ~ Obs 
4/64 - 10/64, 4/65 
4/65, 7/65 


Calc < Obs 
1/65,7/65, 10/65 
10/65 - 7/66 


Thus a, all four 

HoweverX Se« paraph discusses L period of high altitude samplmg a, five latttudes, 
which indicates faster loss from the stratosphere. 

Long-time Observations 

At the conclusion of atmospheric carbon- 14 monitori "g ^f t f970 A Xost ^ears 

ss^sasss5^fe«=sr-«* 

km and up to 36 km at 30N (Telegadas et al„ 1972). 

“P io P r f SSU d al fr Ude i9°53 2 .o 
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0f K arb °, n ; 14 in ^ ho1 ^ air sam P les ’ the AE C developed a molecular sieve collector specific for 
carbon-14 in carbon dioxide. This detector was tested in the AEC high altitude test chamber 
, u a ni der , rea listic stratospheric conduions and tested in the stratosphere with a series of balloon 
launches at 3 IN latitude in 1968 and 1969. Molecular sieve instruments gave carbon- 14 values 

‘ n m a f em , ent . WU ^ n- 7 a n U ?n^ Und ? y oV h L ° le air sam P lin g- A EC conducted an atmospheric research 
program 1 ^'70- 197 1 with 84 balloon launches from six latitudes, 65N, 42N, 3 IN, 9N 

23b, and 34S. Carbon dioxide collected by the molecular sieve absorbers had specific carbon- 14 
two to four times greater than that of carbon dioxide derived from pre-nuclear age standards. In 
multiples of 10 atoms of carbon- 14 per gram of dry air, the natural background of stratospheric 
carbon- 14 was 74 units, and the values observed between 20 and 36 km in 1970-1971 varied 
between 75 and 158 units after the background of 74 units was subtracted (Telegadas, 1971 
Telegadas et al., 1972) This procedure of subtracting 74 units from the total measure of carbon- 
14 is correct only in the troposphere. Jackman et al., 1991, derived the pre-bomb background 
arbon-H using a galactic cosmic ray production function peaking in the stratosphere. This 
study used the magnitude of carbon- 14 produced per year from Lingenfelter and Ramaty (1970) 
but with the altitude distribution correlated with the beryllium-7 production function. Their 
preliminary results suggest that the maximum carbon-14 produced between 20 and 30 km was 76 
o 82 units (Jackman, personal communication). If this is accurate, subtracting a constant 
amoant fo . r P£r' bo ™ b ° r background carbon- 14, at all altitudes and latitudes, would not 
significantly affect the shape of the excess carbon- 14 vertical profiles. Further work should be 
conducted using a production function explicitly for carbon- 14. 

(197 ? discu , ssed whether the Chin ese and French atmospheric nuclear-bomb 
I CS incn *^7-1970 are adequately included in the 1970 observations. During the period of 1953 
o 1959 tb e nuclear _bombs that were exploded in the atmosphere were equivalent in energy to 90 
megatons (MT) of TNT, almost every bomb was less than 10 MT. During 1961-1962 the bombs 
tested were equivalent to 337 MT, including 15 bombs greater than 10 MT, 7 greater than 20 
MT, and one of about 60 MT. The bombs having energy equivalents greater than 10 MT rise 
weU into the ^stratosphere^ After December 1962, atmospheric bomb tests were made as follows: 

4 ? N ’ 3 in 1967 ’ 3 MT in 1968, and in ° ctober 1970; French, 22S, 2.5 and 1.0 MT 
in 1968 and tests in summer 1970. For 3-MT bombs at these latitudes the top of the initial bomb 
cloud is about 22 km. Telegadas compared the atmospheric records of zirconium-95 having 65- 
aay hall-hte with respect to radioactive decomposition and cesium- 137 having a 30-vear half- 
life, and he concluded that the debris of nuclear bomb tests of 1967 to 1970 deposited above 20 

ur 7 i Wa u° ne tW ° orde . rs magnitude less than that deposited below that altitude. In 1970- 
19/1, the maximum mixing ratio of carbon- 14 was observed at altitudes of 24 to 28 km and the 
maximum mixing ratio of zirconium-95 and cesium-137 was between 14 and 18 km in the 
northern hemisphere and between 15 and 19 km in the southern hemisphere. Telegadas 
cautiously concluded that almost all excess carbon- 14 observed above about 21 km in 1970-71 
was from pre-1963 nuclear bomb tests. The results are reproduced in Table 1-7; the calculated 
and observed vertical profiles are given in Figure 1-8; Telegadas 1 contour plot is shown in Figure 
MO The observed vertical profiles are quantitatively different from all the calculated profiles. 
Each observed vertical profile shows a broad maximum carbon- 14 mixing ratio at more-or-less 
Zj km and a distinct decrease with altitude above the maximum; but except for ITALY at 70N 
each calculated profile increases with altitude up to 40 km (Figure 1-8) In the 22- to 28-km 
r^ge NGAR, MPI, GSFC, AER, CALJPL, LLNL, and LLNLND calcula"e carbo„” 4 to be 
much less than the observed values. For example, the observed and calculated values at 26 km 
altitude are listed in Table 1-8. The excess carbon- 14 in the troposphere was 50 units, and the 
ment of a model in predicting long-term stratospheric carbon- 14 is the difference between the 
str at 0 spher,c value and the tropospheric value. These differences range from 18 for the NCAR 
model to 77 for the ITALY model, compared with 90 observed at 26 km. The DUPONT and 
\a/ a cu m °T e ^ S A a w^ W ^ ith the magnitude of the 26-km data to 75% or better. The LLNLND 
WASH, and CAMED models account for better than 50% of the long-term carbon- 14. The 
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NCAR, MPI, GSFC, AER CALJPL, and LLNL models account for 
carbon- 14. The data show far greater persistence of carbon- 14 in 
shown by any model except ITALY and DUPONT. 


40% or less of the long-term 
the 20- to 25-km range than 


The LLNLND two-dimensional model was used to investigate Telegadas s conclusio 
the French and Chinese tests did not enhance the middle stratosphere observed carbon- 
profiles in November 1970. The model was integrated from October 1963 (using Johnston mitia 
conditions) through November 1970 including the effects of the French and 9 hl " es ® tes [f ‘ 
amount of carbon- 14 from these tests was derived by knowing the time, iatttude, bomb cloud 
stabilization height, magnitude of each nuclear test, distribution within the cloud [Bauer 19 ], 

!Tount If carbon-14 produced per MT-2.0 x 1026 atoms/MT [Telegadas, 1971] The 
results of this study are shown in Figure 1-9. The additional carbon-14 produced from the French 
and Chinese tests increased the carbon- 14 mixing ratio in the northern hemisphere ^iow^Okrm 
The LLNLND model with the additional carbon-14 compares better to the ^observed data tin the 
lower stratosphere It does not explain the observed carbon- 14 profile between 25 and 35 kirn 
Including the French and Chinese nuclear tests has a small impact on the vertical profile of 
carbon- 14 in the southern hemisphere. 


Fieure 1-10 is a latitude-altitude contour map of carbon- 14 mixing ratio units during 
November 1970. It is transcribed from Telegadas' (1972), with a different vertical sc ale, and l ^ 
k-icpH nn the data in Table 1-7 The data for the southern hemisphere are sparse, the northern 
hemtohere it taduTfour Ic rtical profiles. In .he 23- to 30 : kra altitude band the data show a 
minimum with respect to latitude above the tropics. This minimum above the tropics can be 

uni!ke 

whhaUitude is unknown to the authors of this chapter. However, see quotauons from Rev, ewers 
A and B in the previous section. 


MODELING STRONTIUM-90 

As with carbon- 14, strontium-90 data were measured by balloonsondes and a very large 
number of aircraft flights. Telegadas (1967) averaged the data over a 3-month period a £ d 
produced altitude-latifude contour plots. Johnston (1989) used these plots to Produce the 
October 1964 initial conditions for this study. Strontium-90 data like ca ™°n * 4 ’ a ^ e 

nronortional to mixing ratios and were modeled in this manner. Typically, four latitu ( , 

31N, 9N, and 34S) were measured from October 1964 to January 1967 with an altitude range 
from the surface to 30 to 36 km. 

Since strontium-90 rapidly coalesces on aerosol particles, and the settling velocity of ae ™ s °J 
particles is proportional to the square of their radius, the removal of strontium-90 is dependent on 
the size of aerosol particles in the atmosphere. In situ measurements of aerosol size from t e 
earlv 1960s are sporadic. Mossop (1964) reported the size distribution of aerosol particles at 20 
km from early 1963 to July 1964. Measurements show that although large size particles with 
diameters equal to 8 qm were collected in early 1963, the median diameters of aeroso l partic les 
gradually reduced to^O.2 qm by July 1964. Volcanic dust from the eruption c >f J^Agung on 
March 17 1963 may have settled to the lower stratosphere by the end of 1964 After 1963, no 
majo^ volcank eruptions were reported until the eruption of Mt. 

this radionuclides study is to simulate the concentration of strontium-90 from October 1964 to 
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die stmosphere^ “* aSSUme tha ' “ ,hiS p ' ri °« *« ™ only background aerosol panicles in 

proX S sTae P rosds 8'obal coverage of the 

aerosol extinction measured by SAGE II can be used m h ? nstrated the multi- wavelength 
Analysis of SAGE II data shows the lUinTofLm^. ae [ oso1 Slze (Yue et al > 1986). 

decreased since the eruption of Ruiz inNovembe? igA f T" the . Stratos P here has gradually 
particles have returned to background levels bv 1Q8Q w ° We assume th e aerosol 

SAGE II data to deduce the glo!£^ a " d 1989 

calculated as a function of altitude and latitude Th? ^1 , tthng velocities were then 

mass mean radii that conserve the total mass of aerosols i n -I" 8 ve * oclt .!f s were derived from the 
decreases with height due to an increase of ambient temperature IhiTmnd 0 ? ' Sl " Ce aer ° So1 size 
-times that the aerosol becomes small enough to acfl^ 

withAe'obw^^ profnefcataloged\y^hnstoifn9^)'* e Fn l * er ' V h l * ^ ' he m0dds are “"Pared 

for strontium-90 dinted as a gas and is F ? ri T h " me P*™* 1 - r “ults are shown 

in this study. As previous work (Johnston et al 1 Q 7 A i^ e ° r seven m °d e is that participated 
90 profiles derived by the Sels c do not reDresenf f h;.K niS ° n ,’ * 989) SUggested ’ *e snontium- 
not incorporated into the scenario. Large demons ^between'The d ! ! J Setding velocit y is 
upper stratosphere are observed (e.g.. Figure I 11 AnriM Q^ w£ data V* the middle to 
used, the agreement between model-derived and nh 9 ? 6 ‘ Wben the settlln g velocity is 
comparable to that between the nuSdels and cafbom 14 d* a dlstnbutl °"s of strontium-90 is 

are shown for the time' thro * h*rT*h 16 ; b0 km) for strontium-90 
stratospheric residence times are not available* ® cto ^r 1966. Observed global 

wero calculated for the four SeVTn hi KSST <198 ? ) local residence times 
1.6 years. Most of the models derived global V3 UeS ranged from 10 to 

settling velocities were included. Without settling v Jirk v esidence tlmes m this range when 
residence times was from 1.8 (AER and GSFC) to 3.9 yearVaTALY)" 86 ° f gl ° bal stratos P heric 

observed strontium-90 ^eVre^entedfn Table I n^fa 8 ’ th ? ra J los of calculated strontium-90 to 
latitudes, and for each of eight times (when Lh J h*! ° f ^ m0de,S ’ for each of four 
the altitude is that of the obsemd ^maximum t WCre reported) ‘ At each latitude, 

64N, 20 km at 30N and 34S, and 24 km at 9N In th^ h 9 ° m,X "V g rat L 10S ’ which m 17 km at 
a systematic decrease of this ratio with time indicates ml^el remn ° f ) C ^ bon_ 1 4 ’ 11 was f °und that 
observed. Since high altitude (above 2oTm)^ £ m SS remova, 1 of ^ acer gas faster than that 
a faster local removal of tracer ^oukT bt due^remti T' y taken at one latitude ’ sucb 

redistribution of tracer within the stratosphere Wkff fr ° 1 m , tbe stra tosphere or it could be 
altitude strontium-90 measuremenrs c f' W th l he availability of four latitudes with high 
transfer of strontium ttMh^upror stratosphere^ beuprc«* e T"" 8 P ? ve,,ts a "y significam 

al.ema.ives is more likely to be correct. The individual nS«toS‘£f ich ° f 

four latitudes CmietlOT f ccrease stron S i V with time at all 

lower stratosphere much faster than ^at obsemed i5 ' | X Cl sweeps strontium-90 out of the 
shows the same pattern. The CALJPL model shows thp r e / S , extrerne extent, the GSFC model 
and to remain about constant at 64N and in the southern h lc/ob ? ratl ° ^ decrease at 30N and 9N 

t a. these three models remove strontium-90 from the stratospterofas^thanS obse^ed"" ^ 
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trends, the average value of the calc/obs ratio is close to one. 
strontium-90 from the stratosphere at the same rate as that observed. 

inte ^s — £ sis 

than the measured values. 

The strontium-90 data test two , thmgs 

-e nai ftL me « thts 

stratospheric-uopospheric exchange is important ^‘STnot apply to that 

raid btm”me P thrne C SS ape" wS whh the stron.ium-90 data, CAMED, WASH and 
problem. Ine three mooeis uidi ■ A 14 data especially the long-term and middle 

~ ,he 

lower and middle stratosphere. 


CONCLUSIONS 

of air motions and aerosol settling. 

average residence time is 4.9 years and the range is 4.3 to 5.4 years. 

lf the 55 ” y n 1MKS 

errors or real differences oetween . , jo removing carbon- 14 from the area in 

Over an 8-year time period, all models but two remove carbon- 14 from 20 to 30 km altitude 
in the stratosphere much more rapidly than observed. 

two-dimensional models 
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Table 1*1. Carbon- 14 stratospheric residence times in years based on linear least squares 
regression analysis. For each period the model derived burden was integrated globally 
between 16 and 60 km. 


Model 


Jan 65-Jul 66 

Oct 63-Jul 66 

Oct 63-Jan 7 1 

Jan 65-Jan 71 

AER 

1.7 

3.5 

2.7 

3.9 

4.5 

CALJPL 

1.7 

3.9 

2.8 

4.2 

4.8 

CAMED 

2.4 

3.4 

3.4 

4.5 

5.0 

DUPONT 

2.7 

4.9 

3.6 

4.9 

5.4 

GISS 

1.3 

3.2 

2.3 

NA 

NA 

GSFC 

1.7 

2.9 

2.4 

3.8 

4.3 

ITALY 

3.2 

4.4 

4.1 

4.9 

5.1 

LLNL 

2.7 

3.5 

3.1 

4.1 

4.7 

LLNLND 

2.0 

3.6 

2.9 

4.6 

5.3 

MPI 

1.3 

3.5 

2.4 

4.3 

5.3 

NCAR 

1.9 

3.0 

2.5 

3.8 

4.4 

WASH 

2.0 

3.7 

3.4 

4.9 

5.4 


Table I-2a. Grouping of global stratospheric residence times (between 16-60 km) for 
October 1963 through October 1964 


Short 

Medium 

Long 

AER 

CAMED-theta 

ITALY 

CALJPL 

LLNL 


GISS 

GSFC 

LLNLND 

MPI 

NCAR 

WASH 

DUPONT 


1.3- 2.0 

2.4 - 2.7 

3.2 


Table I- 2b. Grouping of global stratospheric residence times (between 16-60 km) for 
January 1965 through July 1966 


Short 

Medium 

Long 

GISS 

AER 

DUPONT 

GSFC 

CALJPL 

ITALY 

NCAR 

CAMED-theta 

LLNL 

LLNLND 

MPI 

WASH 


2.9 - 3.2 

3.4 - 3.9 

4.4 - 4.9 
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Table I-2c. Grouping of global stratospheric residence times (between 16-60 km) for 
October 1963 through July 1966 


Short 

Medium 

Long 

AER 

CAMED- theta 

ITALY 

CAUPL 

DUPONT 


GISS 

LLNL 


GSFC 

WASH 


LLNLND 



MPI 



NCAR 



2.3 - 2.9 

3.1 - 3.6 

4.1 


Table I-2d. Grouping of global stratospheric residence times (between 16-60 km) for 
October 1963 through January 1971 


Short 

Medium 

Long 

AER 

CAMED-theta DUPONT 

CAUPL 

LLNLND 

ITALY 

GSFC 

MPI 

WASH 

LLNL 



NCAR 



3.8 -4.2 

4.3 - 4.6 

4.9 


Table I-2e. Grouping of global stratospheric residence times (between 16-60 km) for 
January 1965 through January 1971 


Short 

Medium 

Long 

AER 

CALJPL 

DUPONT 

GSFC 

CAMED-theta ITALY j 

NCAR 

LLNL 

LLNLND 



MPI 



WASH 

4.3 - 4.5 

4.7 - 5.0 

5.1 -5.4 
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_!Z“ 305 0.39 0.49 0.59 0.56 0.66 0.75 0.75 0 75 0.79 0 90 116 0 98 

*2£2E 0.56 0.70 0.75 0.79 0.97 0.87 0.93 1.11 1.Q8 1.13 1.23 1,36 











Table 1-5. Altitude/km of maximum carbon- 14 mixing ratio from January 1964 to July 1966 
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Table 1-6. Change of carbon- 14 mixing ratio with altitude between 25 and 33 km, (du/dz) 
expressed in carbon- 14 mixing ratio units per kilometer, 3 1 north latitude 


DATE 

OBS 

ITALY 

CAUPL 

1/64 

-32 

-32 

-32 

4/64 

-40 

-44 

-18 

7/64 

-40 

-44 

-11 

10/64 

-32 

-22 

-2 

1/65 

-24 

-12 

+2 

4/65 

-11 

-5 

+4 

7/65 

-12 

-6 

+6 

10/65 

-10 

+4 

+7 

1/66 

-6 

+8 

+6 

4/66 

+5 

+9 

+5 

7/66 

-2 

+9 

+5 


Table 1-7. Excess carbon- 14 mixing-ratio units measured as a function 
altitude for November 1970 [Telegadas et al., 1972] 


of latitude and 


65N 42N 30N 9N 34S 


z 

C 

z 

C 

Z 

31.0 

120 

36.0 

125 

36.3 

30.0 

130 

31.5 

130 

35.5 

27.0 

140 

30.8 

134 

34.0 

26.9 

141 

27.4 

138 

33.0 

24.0 

139 

27.2 

140 

32.8 

22.3 

140 

24.0 

143 

32.5 

20.3 

130 

23.3 

140 

31.2 



21.3 

130 

27.6 



20.9 

125 

27.2 



19.6 

120 

26.5 


24.4 

24.3 
23.9 
22.2 
21.2 
21.0 

20.3 
19.6 


C 

z 

C 

Z 

C 

96 

31.5 

106 

32.3 

120 

110 

30.6 

no 

27.3 

120 

120 

27.2 

in 

27.2 

130 

130 

23.7 

120 

27.0 

134 

128 

23.0 

110 

24.2 

135 

116 

22.0 

100 

23.9 

130 

130 

21.4 

90 

21.7 

120 

142 

21.0 

80 

21.1 

128 

131 

20.8 

77 

20.9 

119 

140 



20.0 

120 


139 

140 
142 
130 
128 
135 
120 
110 
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Table 1-9. Settling velocities (x l<Hcm s' 1 ) of aerosol particles used in the strontium-90 simulation 
Southern Hemisphere Latitude Zones in Degrees 
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2.5 282.5 





globally from 16 to 60 km. 


AER 
CAUPL 
CAMED-theta 
GSFC 
ITALY 
LLNL 
WASH 


0.9 

1.4 
1.6 
1.1 
2.1 

1.5 

1.6 


1.8 

NA 

2.9 

1.8 

3.9 
2.3 
2.3 
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Table 1-11 The ratio 90 Sr (calcar (obs) at four latitudes and at the altitude of maximum 

observed yu Sr at each latitude for calculations that included settling velocities of stratospheric 
Junge particles 
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FIGURE CAPTIONS 


Figure 1-1 Carbon-14 altitude-latitude initial conditions for October 1963: a) Telegadas (1971) 
constructed contours for October 1963 (which is an average of September .October and 
November data); b) Johnston (1989) initial conditions for October 1963 (used m this study), and 
c) Kinnison (1989) model-derived October 1963 data. Contours are in units of 105 atoms of 
carbon- 14 per gram of dry air as obtained from a total air sampler. These units are proportional 

to mixing ratio. 

Figure 1-2, Individual measurements at 31N of ^ton-14 over a three-month time penod 
centered on October 1963, January 1964, October 1964, and October 1965. The solid line for 
each time period is the profile derived by Johnston (1989) from Telegadas contours. 


Figure 1-3. Comparison of model-derived carbon- 14 and observed data at 70N between January 
1964 and July 1966. Profiles are in units of 105 atoms of carbon- 14 per gram of dry air as 
obtained from a total air sampler. These units are proportional to mixing ratio. 


Figure 1-4. Comparison of model-derived carbon- 14 and observed data at 31N between January 
1964 and July 1966. Profiles are in units of 105 atoms of carbon-14 per gram of dry air as 
obtained from a total air sampler. These units are proportional to mixing ratio. 

Figure 1-5. Comparison of model-derived carbon- 14 and observed data at 9N between January 
1964 and July 1966. Profiles are in units of 105 atoms of carbon- 14 per gram of dry air as 
obtained from a total air sampler. These units are proportional to mixing ratio. 


Figure 1-6. Comparison of model-derived carbon-14 and observed data at 42S between January 
1964 and July 1966. Profiles are in units of 105 atoms of carbon- 14 per gram of dry air as 
obtained from a total air sampler. These units are proportional to mixing ratio. 

Figure 1-7. Model-derived carbon-14 altitude-latitude contour plots for October 1964 and 
October 1965. Contours are in units of 105 atoms of carbon- 14 per gram of dry air as obtained 
from a total air sampler. These units are proportional to mixing ratio. 

Figure 1-8. Comparison of model-derived carbon-14 and observed data at 70N, 42N, 30N, 9N, 
and 34S for November 1970. Profiles are in units of 105 atoms of carbon-14 per gram of dry air 
as obtained from a total air sampler. These units are proportional to mixing ratio. 


Figure 1-9. Comparison of LLNLND model-derived carbon- 14 with and without additional 
carbon- 14 from the French and Chinese nuclear tests (November 1970). Profiles are in units of 
105 atoms of carbon- 14 per gram of dry air as obtained from a total air sampler. These units are 
proportional to mixing ratio. 

Figure I- 10. Contour plot of mixing ratio of carbon-14 in November 1970. Equal distance 
along the horizontal scale corresponds to equal surface area on the three-dimensional globe, and 
equal distance along the vertical scale corresponds approximately to equal mass of air. ihe 
dashed line represents a standard tropopause. 

Figure 1-11. Comparison of model-derived strontium-90 and observed data at 64N between 
April 1965 and July 1966. For each time period, results are shown for strontium-90 treated as a 
gas and as an aerosol particle. Profiles are in units of disintegrations per minute per 1000 cubic 
feet (28 m3) of standard air. These units are proportional to mixing ratio. 
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Figure 1-12. Comparison of model-derived strontium-90 and observed data at 31N between 
January 1965 and January 1967. For each time period, results are shown for strontium-90 treated 
as a gas and as an aerosol particle. Profiles are in units of disintegrations per minute per 1000 
cubic feet (28 m3) of standard air. These units are proportional to mixing ratio. 

Figure 1-13. Comparison of model-derived strontium-90 and observed data at 9N between 
January 1965 and January 1967. For each time period, results are shown for strontium-90 treated 
as a gas and as an aerosol particle. Profiles are in units of disintegrations per minute per 1000 
cubic feet (28 m3) of standard air. These units are proportional to mixing ratio. 

Figure 1-14. Comparison of model-derived strontium-90 and observed data at 30S for January 
1965 and 34S for April 1965 through January 1967. For each time period, results are shown for 
strontium-90 treated as a gas and as an aerosol particle. Profiles are in units of disintegrations 
per minute per 1000 cubic feet (28 m3) of standard air. These units are proportional to mixing 
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Carbon-14 Mixing Ratio Units Figure i-5 (cont.) Carbon-14 Mixing Ratio Units 
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Carbon-14 Mixing Ratio Units Figure i-5 (cont.) Carbon-14 Mixing Ratio Units 
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Carbon-14 Mixing Ratio Units Figure i-6 (cont.) Carbon-14 Mixing Ratio Units 
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Carbon-1 4 Mixing Ratio Units Figure i-6 (com.) Carbon-1 4 Mixing Ratio Units 
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Carbon-14 Mixing Ratio Units Figure i-6 (cont.) Carbon-14 Mixing Ratio Units 
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Figure 1-6 (cont.) 
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LATITUDE Figure 1-7 (cont.) LATITUDE 
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Carbon-14 Mixing Ratio Units Figure i-8 Carbon-14 Mixing Ratio Units 
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Carbon-14 Mixing Ratio Units Figure i-8 (cont.) Carbon-14 Mixing Ratio Units 
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Figure 1-8 (cont.) 
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Carbon-14 Mixing Ratio Units Figure i-9 Carbon-14 Mixing Ratio Units 
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Strontium-90 Mixing Ratio Units Figure 1-12 (cont.) Strontium-90 Mixing Ratio Units 
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Strontium-90 Mixing Ratio Units Figure 1-12 (cont.) Strontium-90 Mixing Ratio Units 
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INTRODUCTION 

transnoir P Th^p ae nf ? S i fr ° m the eru , ptio " of a volcano are a unique indicator of stratospheric 
transport These particles may evolve from the S0 2 injected by the volcano. They are 

globally, and eventually are removed by both sedimentation and large-scale transport 

the emotion C I 3! a?. & t0 ,?. few . years ’ de P endin g on the magnitude and location of 

^ eru P tlor V Like the exotic radionuclides, they provide a direct measure of the global spread of 
a species injected in the lower stratosphere. e p or 

The eruption on 13 November 1985 of the volcano Nevado del Ruiz (5N, 75 W) in Columbia 
was observed by the Stratospheric Aerosol and Gas Experiment II (SAGE II) satellite 
experiment. This was the strongest volcanic eruption observed by SAGE II before the erumion 
o inatubo on 15 June 1991 . The optical depth of aerosol particles associated with Ruiz reached 
maximum in Febniary 1986 and about 2 years later decayed into the background levels of 
e osols left by El Chichon (Yue et al., 1991). Initial distributions (z* by latitude) of particle 
number concentrations for 1 February 1986 were derived from the SAGE II data set by assuming 
a lognormal si ze distribution (Yue et al., 1986). In addition, the settling velocities were also 

nSTS by hC S1ZC dl f nbut ‘ ons inferred from ^e SAGE II data set (Kasten, 1968). Each model 
mnH . the sam *r ini “ al conditions and settling velocities in their model simulation and the 
modeling resultsof aerosol number concentrations measured on the 15th of each month were 
intercompared. The capability of simulating the transport of a tracer by two-dimensional models 
was studied by comparing modeling and measurement results. It should be noted that because of 
the possible presence of clouds around the tropopause and much larger uncertainties of aerosol 
extmcuon at shorter wavelengths in the troposphere, the size distributions of aerosol particles in 
TH P ^ POSP i here t and - ar ° Und ? e tr °P°P ause derived from the SAGE II data set are not reliable 
urme ri ai n ri»»e CX j 1 n C at shortei i wavelengths at altitudes higher than 30 km also have large 

aitS k ' If i addition aerosol particles become small enough to act like a gas (v se t = 0) at 
altitudes above 30 km. The comparison among modeling results discussed in this paper is 
limited to altitudes between 10 and 30 km, and the comparison between modeS and 
experimental results is limited to altitudes between the tropopause and 30 km. g 

To avoid the complex level of model development required on most two-dimensional models 
for accurate aerosol microphysical simulations aerosol particles are treated as one constituent 

re/aHteJvflro 126 ' Such treatm ^ nt ls understandably oversimplified: stratospheric aerosols have a 
relatively large range of sizes that continually change due to microphysical processes including 
crnguktion, condensation, evaporation, and sedimentation. However, the change of aerosol size 

wafhemmff 011 d T 3 ' d " amatic since this was a moderate eruption and the simulation 

^nverterto sulfate” 10 " HS aft£r thC erUptl ° n when most of the excess S °2 already had been 
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INTER-MODEL AND MODEL-EXPERIMENT COMPARISONS 


Figure J-l shows altitude-latitude contour plots of aerosol number concentration, based on 
SAGE II observations, used by all models as initial conditions. Model simulation started on 1 
February 1986. The initial concentrations of aerosol particles at altitudes below 8 km are all set 
to 2 per cm^. The peak of Ruiz aerosols is at 20 km and over the equator. The injection occurred 
during the easterly phase of the quasi-biennial oscillation (QBO) when lofting and equatorial 
confinement are expected (Trepte and Hitchman, 1992). 

The initial aerosol number concentrations and settling velocities are listed in Tables J-l and 
J-2, respectively. The settling velocities were derived from the mass mean radii that conserve the 
total mass of aerosols in a size distribution. 


Seven groups have submitted model results for this experiment. The contour plots of aerosol 
number concentration for the months of April and August 1986 are shown in Figures J-2 and J-3, 
respectively The intervals between contours in these plots are 1. It can be seen that among 
seven two-dimensional models, AER and CALJPL have the most rapid transport out of the 
tropics In April, while most models show a maximum of about 5 at 20 km over the equator, the 
maximum is only about 3 in the CALJPL and AER models. The maximum value derived from 
the SAGE II data set is about 8, which is higher than the values obtained by all model 
simulations. In addition, the location of peak concentration has moved up to 22 km and remains 
there for the rest of the year. By August, most of the Ruiz debris has been transported to higher 
latitudes in the AER and CALJPL models. For the other five models, the maximum 
concentration of Ruiz aerosols is still at about 20 or 22 km, and is less than observed but still in 
reasonably good agreement with the experimental result. It is of interest to note that the 
maximum aerosol concentrations for the GSFC and LLNL models are no longer over the equator 
- they have shifted to higher latitudes in the southern hemisphere. 

Profile comparisons over the equator are shown in Figures J-4a, b, c, and d, for the months of 
March, June, September, and December 1986, respectively. For the month of March 1986, the 
differences of aerosol concentrations among models and experimental results are very small lor 
altitudes above 25 km. The small coefficients k zz above 25 km used in the models are in good 
agreement with the weak lofting observed by the experiment. However, in regions around the 
maximum concentration there are large differences between model and observational resu ts. 
Among model results, the CALJPL and AER have the lowest peak concentration values as 
shown in the contour plots. The CALJPL model normally runs with the dynamics specified as a 
streamfunction. These runs were made in a nonstandard mode with dynamics specified as 
velocity fields, and the results may contain a small error. Three months kter better agreement is 
seen between SAGE values and the model results of CAMED, GSFC, ITALY - 1 
WASH, even though below 19 km the CAMED results are higher than observed an d the WASH 
results are lower. Since in the previous months the peaks of profiles of AER and CALJPL are 
lower than the peaks of profiles of other model results, it is not surprising that they remain lower 
for subsequent months. By September WASH results are in good agreement with the 
observations, but the CAMED results below 19 km are much higher. By December 1986 the 
broadening of the peak is shown in all results. The comparison is excellent insofar as the setthng 
velocity is needed to explain the first-order effects on the aerosol. After June the GSFC, ITALY, 
LLNL, and WASH models do reasonably well, indicating possibly a good averaged circulation. 
The CAMED model looks good but is centered too low; perhaps there is a problem with the 
effective tropopause. The AER and CALJPL models diverge significantly over the long term. 
All of the models have difficulty preserving the peak concentrations; they seem to wash out the 
detailed structure. 
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Profile comparisons at 35S and 35N for June 1986 are shown in Figures J-4e and J-4f, 
respectively. Although there are discrepancies about the location and value of the peak among 
different model profiles, they all demonstrate the influence of Ruiz aerosols on the number 
concentrations at altitudes below 22 km. It is of interest to note that the influence of Ruiz 
aerosols does not appear in SAGE II profiles above 20 km at 35N, since the main Ruiz layer was 
confined in tropical regions between 30S and 30N. 

The temporal variations of aerosol number concentration at different locations are shown in 
Figures J-5a, b, c, and d. Figure J-5a shows that near the injection point, in the first few months 
of simulation, the decay rates of all models are faster than the measurement results. However, 
after July the decay rates are all in rather good agreement. The differences in concentrations for 
different models in the later months of ihe year are mainly due to the differences resulting from 
transport in the early months of the simulation. It is of interest to note that already by mid- 
February, the concentration of aerosol particles in the CALJPL model is much lower than that in 
other models. This may be a result of stronger diffusivity used in the model. However, it is 
more likely a basic problem with the initialization to the grid points of each model. In general, 
there are large differences in decay rates among model and observational results in the first few 
months of simulation. However, the decay rates after July 1986 at different altitudes at the 
equator shown by different models and the SAGE II data set are in reasonably good agreement. 
At low latitudes, the agreement between experimental and modeling results is still quite good; 
however, transport to higher latitudes in the models is much faster than observed (Figure J-5d). 

A quantitative assessment of the transport in each model can be conducted by examining the 
model-derived maximum concentrations at several latitudes in the tropical region listed in Table 
J-3. The AER and CALJPL models have more rapid mixing out of the tropics. By May 1986, 
AER results show that the peak concentration of Ruiz aerosols has disappeared from the equator 
and has been transported to higher latitudes in both hemispheres. The same is true for the 
CALJPL model in Junel986. In both models the peak number concentrations have been reduced 
to about 2 by June. On the other hand CAMED, ITALY, and WASH have weak transport in 
comparison with other models. They preserved the peak aerosol concentration in the tropical 
region almost until the end of the year. The transport in the GSFC and LLNL models is in- 
between, since by June the aerosol concentrations over the equator were relatively high, but the 
peak gradually shifted to higher latitudes in the southern hemisphere in subsequent months. 

In Table J-4, the SAGE II and model-derived aerosol concentrations at 0 degree are shown 
for altitudes from 16 to 24 km. Since aerosol particles over the equator are rapidly transported to 
higher latitudes in the AER and CALJPL models, the effect of upwelling on aerosol particles is 
weak. On the other hand the effects of upwelling can be seen in GSFC and ITALY models, since 
the peaks of aerosol concentration have been gradually shifted to 22 km. For the WASH model, 
the concentrations in the altitude range 22 and 24 km are high in comparison with the very low 
concentrations below 20 km. Additionally, the peak of aerosol concentration has shifted to 24 
km by the end of 1986. It is of interest to see that the CAMED model not only has weak 
transport but also weak upwelling. By the later months of the year, relatively high 
concentrations of aerosol particles have remained at all altitudes below 22 km in the CAMED 
results, while all other models show very low concentrations due to either rapid advection or 
diffusion. 

It should be noted that all model results have been obtained with parameterized washout 
included. We have tried to compare model results with relatively large initial concentrations of 
aerosol particles in the troposphere and no washout processes. It has been found that some of the 
aerosol number concentrations in the lower stratosphere are unrealistically high. Clearly we need 
more accurate microphysical models for removal of aerosols near the tropopause if we are to 
include tropospheric aerosols in these simulations. 
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CONCLUSIONS 


Comparison between models and observations has shown that all models do poorly in the 
first few months of simulation when there are large aerosol gradients. When gradients have 
become smoother a few months later, the models perform better. Intercomparison of the models 
shows that the AER and CALJPL models have the strongest transport out of the tropics, losing 
the peak concentrations in just a few months. The peaks of aerosol concentrations remained over 
the equator in the CAMED, ITALY, and WASH models, whereas the GSFC and LLNL models 
have stronger transport in the southern hemisphere that developed high-latitude maximum 
concentration through transport. The GSFC, ITALY, and WASH models seem to have strong 
upwelling, whereas the CAMED model has less upwelling. Alternatively, there may be 
problems with the interpolation of the initial conditions in some of the models. 

A given eruption occurs during a specific phase of the QBO. During easterly shear, aerosols 
will be lofted and confined in the tropics, whereas during westerly shear, aerosols may spread out 
more easily, being more rapidly transported poleward and downward a few kilometers above the 
tropopause. SAGE II results show that Ruiz aerosols are confined in the tropics, but this behavior 
is not being taken into account in current models. The comparison between models and the 
observed evolution of volcanic debris suggests that the influence of the QBO on the mean 
tropical circulation may need to be included in future model development. This circulation may 
be specified readily from previous theoretical and modeling studies of the QBO. 

The good agreement of decay rates obtained by modeling and measurement results in the 
later months of simulation suggests that modeling can be a useful tool to estimate the global 
decay of volcanic debris or engine exhaust from supersonic aircraft in the stratosphere. On the 
other hand, the accuracy of simulating the advective and diffusive processes in the models can be 
improved by comparing modeling and experimental results when the opportunity occurs. Better 
agreement is expected if we can have a better dynamical model for the exchange and transport 
near the tropopause, if it can be modified to have better vertical and latitudinal resolutions, and if 
it can include some microphysical processes that are taking place in the volcanic aerosols. 
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Table J-l: Initial Concentrations (number/cm 3 ) at Seven Latitudes 


Z* (km) 

30°S 

20°S 

10°S 

0° 

10°N 

20°N 

30°N 

28 

1.5 

1.3 

0.9 

1.1 

0.8 

0.9 

1.2 

26 

1.4 

1.4 

1.2 

2.5 

1.1 

0.9 

1.0 

24 

1.0 

1.3 

1.5 

4.2 

2.4 

1.1 

1.2 

22 

1.5 

2.2 

4.1 

8.6 

4.1 

1.5 

1.6 

20 

2.6 

4.6 

7.1 

10.8 

6.8 

4.1 

2.4 

18 

4.7 

4.9 

6.3 

5.9 

7.4 

9.1 

4.9 

16 

5.1 

4.6 

3.8 

3.3 

3.0 

6.3 

6.2 

14 

3.2 

3.0 

4.3 

4.0 

4.0 

2.6 

4.3 


Table J-2: Initial setting velocities (cm/sec) at Seven Latitudes 


7 * (km) 

30°S 

20°S 

10°S 

0° 

10°N 

20°N 

30°N 

28 

2.2e-2 

2.5e-2 

3.4e-2 

3.6e-2 

3.4e-2 

2.7e-2 

2.3e-2 

26 

1.9e-2 

2.1e-2 

2.7e-2 

2.8e-2 

2.7e-2 

2.2e-2 

2.1e-2 

24 

1.8e-2 

1.8e-2 

2.1e-2 

2.1e-2 

2.1e-2 

1.8e-2 

1.7e-2 

22 

1.5e-2 

1.5e-2 

1.6e-2 

1.6e-2 

1.6e-2 

1.5e-2 

1.4e-2 

20 

l.le-2 

1.0e-2 

1.2e-2 

1.2e-2 

l.le-2 

l.le-2 

l.le-2 

18 

7.2e-3 

7.3e-3 

7.3e-3 

8.3e-3 

7.4e-3 

6.1e-3 

8.0e-3 

16 

5.0e-3 

5.2e-3 

5.4e-3 

6.6e-3 

5.8e-3 

4.6e-3 

5.7e-3 

14 

4.5e-3 

4.3e-3 

4.2e-3 

5.4e-3 

4.6e-3 

4.5e-3 

5.0e-3 
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-3: Model-Derived Maximum Concentrations (number / cm3) at Seven Latitudes 

” AER 1 

lhn°q l?0°sl 10°S I 0° I 10°N| 20°N| 30°N I 


Month 30°S 20°S 

February 3.1 3.7 

March 2.4 3.1 

April 2.3 2.5 

May 2.3 2.2 

June 2.2 2.1 

July 2.2 2.0 

August 2.1 1.9 

September 1.9 1.7 

October 1.8 1.6 

November 1.6 1-3 

December 1.5 1.1 


10°S 0° 

5.4 6.9 

3.9 4.5 

2.8 3.1 

2.3 2.4 
2.0 2.0 

1.8 1.8 

1.6 L6~ 

1.4 1.4 

1.2 1.2 
1.0 _ 1 £_ 
0.9 0.9 


20°N 

4.0 

3.9 

3.5 

3.0 
2.5~ 

2.1 

1.9 
1.7 

1.6 
1.6 
1.5 


30°N 

4.2 

3.7 

3.6 

3.2 
2.9 
2.5 

2.3 
2.1 
2.0 

1.8 

1.7 


CALJPL 


Month 


February 


March 


April 


June 


July 


August 


September 


October 


November 


3.2 3.7 


2.8 3.0 


2.5 2.6 


10°S 

0° 

10°N 

20°N 

KB 

4.4 

4.2 

KB 

3.6 


3.5 

3.0 


3.2 

3.0 

2.6 

2.6 

2.6 

2.5 

2.3 

2.2 

2.2 

2.1 

2.0 

1.9 

ia 

1.8 


■H 

1.7 

1.6 

1.6 


.6 1.6 


.5 1.4 


.3 


1.1 1.1 


-6 




































































Table J-3 (cont.) 


CAMED 


Month 

30°S 


0° 


30°N 

February 

4.3 

4.3 

5.6 

PI 

7.0 

5.3 

5.0 

March 

3.9 

4.2 

5.4 

6.6 

6.1 

5.0 

4.7 

April 


5.7 

5.5 

4.7 

4.4 

May 

3.6 

4.0 

4.8 

5.2 

5.0 

4.4 

4.1 

June 

3.7 

3.9 

4.5 

4.7 

4.6 

4.1 

! 3.9 

July 

3.7 

3.9 

4.2 

4.4 

4.3 

3.9 

3.7 

August 

3.6 

3.8 

4.0 

EB 

4.0 

3.7 

3.6 

September 

3.5 

3.6 

3.6 

3.8 

3.8 

3.6 

3.5 

October 

3.3 

3.4 

3.5 

3.6 

3.6 

3.5 

3.5 

November 

3.1 

3.2 

3.3 

3.4 

3.4 

3.3 

3.5 

December 

2.9 

3.0 

3.1 

3.2 

3.3 

3.3 

3.5 


GSFC 


Month 

| 30°S 20°S 

10°S 

0° 

10°N 

20°N 

30°N 

February 

IBS 

IKS 

6.7 8.7 

6.6 

5.5 

4.6 

March 

3.6 

4.1 

5.7 

6.5 

5.1 

4.3 

4.0 

April 

3.3 

3.9 

5.0 

5.2 

4.3 

3.6 

3.6 

May 


3.7 

4.3 

4.3 

3.7 

3.3 

3.3 

June 

pi 

m 

3.8 

3.8 

3.3 

3.0 

2.4 

July 

3.0 

3.2 

3.5 

3.4 

2.9 

2.6 

2.7 

August 

3.0 

3.1 

3.2 

3.0 

2.6 

2.3 

2.4 

September 

2.8 

2.8 

2.9 

2.6 

2.3 

El 

2.2 

October 

2.5 

2.5 

2.6 

2.4 

IBI 

2.0 

2.1 

November 

2.3 

2.3 

2.3 

2.1 

1.9 

1.9 

2.0 

December 

2.0 

IBI 

2.0 

m 

1.8 

1.9 

1.9 













































































































































































Table J-3 (cont.) 
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Table J-3 (cont.) 


WASH 

Month 

m 


10°S 

0° 

10°N 

20°N 

30°N 

February 

4.2 

4.0 

4.6 


5.8 

4.6 

4.8 

March 

3.8 

3.6 

4.7 

6.3 

4.6 

3.4 

4.4 

mm 

3.4 

3.3 

4.6 

5.6 

4.7 

3.0 

3.5 

| 

3.2 

3.0 

4.4 

5.0 

4.0 

2.8 

2.9 


R$li 




Blgl 

2.8 

2.7 

mum 

n 




IBI 


B19 

mi 



ms 

m 

■II 

m 

II 



2.9 

3.5 

3.6 

3.2 

2.9 

2.9 

October 

2.7 

2.8 

3.3 

3.3 


2.8 

2.7 

November 

2.4 

2.6 

2.9 

3.0 

2.9 

2.6 

2.6 

December 

IBI 

tsa 

2.5 

2.6 

2.8 

2.5 

2.1 
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Table J-4: SAGE II and Model-Derived Concentrations ( number/ cm 3 ) at 0 C 


SAGE II 


18 km 


20 km 

22 km 

7.9 

9.4 

6.5 

7.9 

5.4 

6.4 

4.2 

4.8 

3.8 

4.5 

3.7 

4.4 

3.0 

3.7 

2.7 

3.0 

2.5 

2.5 

2.3 

2.3 


Month 


March 


April 


May 


June 


July 


August 


September 


October 


November 


December 



16 km 


2.8 


2.4 


2.3 


2.2 


2.0 


2.0 


.9 


.9 


1.8 


24 km 


.6 



Month 


February 


March 


April 


May 


June 


July 


August 


September 


October 


November 


December 


16 km 


3.0 


1.9 


.3 


0.9 


0.7 


0.6 


0.5 


0.4 


0.4 


0.3 


0.3 


18 km 


5.9 


3. 


2.4 


1.8 


1.5 


.3 


1.1 


0.9 


0.8 


0.7 


0.7 


20 km 

22 km 

6.9 

5.7 

4.5 

4.0 

3.1 

3.0 

2.3 

2.4 

1.9 

2.0 

1.7 

1.8 

1.4 

1.6 

1.2 

1.4 

1.1 

1.2 

1.0 

1.0 

0.9 

0.9 


24 km 


4.0 


3.0 


2.4 


2.0 
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Table J-4 (cont.) 


CALJPL 


Month 


February 


March 


April 


May 


June 


July 



October 


November 


December 


16 km 


1.6 



18 km 


2.8 


.9 


1.3 


1.0 


0.9 


0.8 


0.7 


0.6 


0.6 


0.5 


0.4 


20 km 


4.4 


3.4 


2 . 


2.2 


22 km 


4.4 


3.7 


3.2 


2.6 




24 km 


3.3 


2.9 


2.5 


2.2 



CAMED 


Month 


March 


A 


June 


September 


October 


November 


December 


16 km 


4.6 


4.2 


3.6 


3.3 



18 km 


6.0 


5.5 


5.0 


4.6 


4.2 


4.0 


3.8 


3.5 


3.3 


20 km 

22 km 

7.9 

6.2 

6.6 

5.0 

5.7 

4.4 

5.2 

3.9 

4.7 

3.6 

4.4 

3.3 

4.1 

3.1 

3.8 

2.9 

3.6 

2.7 

3.4 

2.5 

3.2 

2.4 


24 km 


3.5 


3.0 


2 . 


2.5 


2.3 


2.2 


2.0 


1.9 


1.8 



J-ll 




























































































































Table J-4 (cont.) 




J- 12 












































































































































Table J-4 (cont.) 


October 


November 


December 


LLNL 


20 km 


6.0 


5.0 


4.3 


3.8 


3.5 



1.1 

2.6 

3.0 

2.7 

0.9 

2.3 

2.7 

2.5 

0.7 

msm 

2.4 

2.3 

0.6 

1.7 

2.2 

2.1 



WASH 


Month 


March 


April 



une 


July 


August 


October 


November 


16 km 

18 km 

20 km 

4.2 

5.2 

7.2 

2.5 

2.7 

5.3 

1.4 

1.2 

3.6 

1.0 

0.7 

2.7 

0.8 

0.8 

2.4 

0.8 

0.9 

2.4 

0.8 

0.9 

2.4 

0.5 

0.8 

2.3 

0.4 

0.6 

1.9 

0.4 

0.6 

1.6 

0.4 

0.5 

1.2 


22 km 


6.9 


6.3 


5.6 


5.0 


4.3 


4.0 


3.8 


3.6 


3.3 


3.0 


24 km 


4.3 
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FIGURE CAPTIONS 

Figure J-l. Initial concentrations (number/cm3) for model simulation. 

Figure J-2. SAGE II and model-derived concentration (number/cm3) contour plots for April 
1986. 

Figure J-3. SAGE II and model-derived concentration (number/cm3) contour plots for August 
1986. 

Figure J-4. Comparison of profiles of aerosol number concentration. 

(a) 0°, March 1986; 

(b) 0°, June 1986; 

(c) 0°, September 1986; 

(d) 0°, December 1986; 

(e) 35S, June 1986; 

(f) 35N, June 1986. 

Figure J-5. Comparison of temporal variations of aerosol number concentration. 

(a) 20 km, 0°; 

(b) 26 km, 0°; 

(c) 20 km, 15S; 

(d) 20 km, 25N. 
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K: Photodissociation Rates 



Overview 

Sections K, L, and M: The ATMOS Profiles 


Sections K, L, and M have a common theme in that we have adopted a mean atmospheric 
profile based on the ATMOS solar occultation measurements made dunng a Space Shuttle tlignt in 
May 1985 Section K examines photolysis rates (i.e., J-values), and section L compares the 
model chemistry of rates and radicals. These sections are basically model intercompansons; they 
address critical issues that remain unresolved since the last model-model mtercomparison Section 
M is a true model and measurement comparison: the models were meant to simulate the sunset 
measurements made at latitudes near 30N. The ATMOS 30N atmosphere defined in the table 
below is a mixture of ATMOS observations, some model interpolation, some climatologies and 
tracer correlations (see section H), and a bit of guess work. It should “ S™Tof 

framework with which to compare the basic photochemistry of the models independent ot 

transport. 


TABLE KLM-1. The ATMOS 30N Atmospheric Profile 


z* 

P 

T 

O 3 

NOy 

(km) 

(mbar) 

(K) 

(ppm) 

(ppb) 

0 

1000 

294 

0.05 


2 

750 

284 

0.05 


4 

562 

272 

0.05 


6 

422 

259 

0.05 


8 

316 

245 

0.05 


10 

237 

231 

0.10 


12 

178 

220 

0.20 


14 

133 

204.1 

0.34 

0.80 

16 

100 

205.6 

0.72 

1.09 

18 

75 

207.9 

1.38 

2.16 

20 

56.2 

210.8 

2.34 

4.24 

22 

42.2 

214.1 

3.53 

6.48 

24 

31.6 

217.8 

4.88 

8.59 

26 

23.7 

221.9 

6.18 

10.83 

28 

17.8 

226.3 

7.15 

12.89 

30 

13.3 

230.7 

7.76 

15.15 

32 

10.0 

234.3 

8.07 

15.93 

34 

7.50 

237.8 

8.08 

15.40 

36 

5.62 

243.1 

7.71 

16.06 

38 

4.22 

249.6 

6.95 

16.20 

40 

3.16 

255.6 

5.87 

16.71 

42 

2.37 

261.3 

4.77 

16.76 

44 

1.78 

267.5 

4.04 

16.91 

46 

1.33 

272.1 

3.53 

15.35 

48 

1.00 

273.1 

3.06 

13.14 

50 

0.750 

270.7 

2.60 

11.11 

52 

0.562 

265.7 

2.13 

8.87 

54 

0.422 

259.6 

1.65 


56 

0.316 

253.3 

1.22 


58 

0.237 

247.2 

0.95 


60 

0.178 

241.0 

0.77 


62 

0.133 

234.7 

0.63 


64 

0.100 

228.7 

0.50 


66 

0.075 

223.4 

0.37 


68 

0.056 

218.8 

0.24 



CH 4 H 2 O Cl y # Aerosol Area* 

(ppm) (ppm) (ppb) ( 10 ' 8 /cm) 


Use latitude = +30° 

Use solar declination = +15° 

Use Br y = 15 ppt everywhere 


1.61 

6.01 

0.20 

1.0 

1.57 

4.14 

0.30 

1.0 

1.50 

3.39 

0.60 

1.0 

1.39 

3.84 

1.20 

0.8 

1.26 

4.62 

1.70 

0.6 

1.12 

4.90 

2.15 

0.4 

0.98 

4.88 

2.45 

0.3 

0.86 

5.07 

2.60 

0.2 

0.77 

5.41 

2.75 

0.1 

0.74 

5.57 

2.80 

0.05 

0.76 

5.53 

2.75 

0 

0.80 

5.47 

2.80 

. 

0.82 

5.48 

2.80 

. 

0.78 

5.60 

2.80 

. 

0.70 

5.78 

2.80 

. 

0.60 

5.99 

2.80 

. 

0.52 

6.24 

2.80 


0.45 

6.50 

2.80 


0.39 

6.63 

2.80 


0.33 

6.60 

2.80 



ft Uses empirical NOy-Cly from ER-2 data 

$ Use k = 5200 x AA x G (sec 4 ) 

G = 0.1 for N2O5 + aerosol 
G = 0.006 exp[-0.15(T-200)] 
for CIONO2 + aerosol 
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K: Model-Model Intercomparison: Photodissociation Rates 


Richard S. Eckman 
NASA-Langley Research Center 

Gail P. Anderson 

Air Force Geophysical Laboratory 
Yuk L. Yung 

California Institute of Technology 


INTRODUCTION AND EXPERIMENT DESCRIPTION 


A disturbing conclusion from the 1988 intercomparison of two-dimensional models 
(Jackman et al., 1989) was that substantial unresolved differences existed, of the order of a factor 
of two, among the photodissociation rates calculated for a specified atmosphere. Not all 
modelers use the same calculation methods or parameterizations, particularly for absorption in 
the O 2 Schumann-Runge band region. Nonetheless, a degree of consistency is essential for an 
intercomparison of other model-calculated chemical quantities to be meaningful. 

For the 1992 Models and Measurements Workshop, further constraints were placed on the 
specification of the model atmosphere. A standard atmosphere based on Atmospheric Trace 
Molecule Spectroscopy (ATMOS) occultation profiles at 30N was used to specify the 
temperature and ozone vertical profiles. The solar zenith angle was fixed at 15 degrees, which is 
the value at local noon for this latitude on May 1 . Four cases were considered: 

Case 1. A full calculation including Rayleigh-phase scattering and surface albedo of 0.3, 
and J-values reported at all solar wavelengths. 

Case 2. Scattering and albedo included, but J-values reported only below 200 nm in an 
effort to isolate the contribution of the O 2 Schumann-Runge band absorption . 

Case 3. No scattering and no surface albedo; J-values reported at all wavelengths. 

Case 4. As above, but J-values reported at wavelengths less than 200 nm only . 

Noontime values were reported in all four cases with a 24-hour average also calculated for 
case (1). Many models use diurnal averages of the photodissociation rates, and it is only the 
average that is important in stratospheric photolysis of O 2 and CFCs. Thus, it was felt that this 
final case would be illustrative of how models perform this average. For a comparison of diurnal 
averaging of species and rates see section L. 

The photodissociation rate (or J-value) is defined as: 

}{z,X)= integral tx(A) F( X,z,X) dX 

where a is the cross section of the molecule, and F(A,z,Z) is the sum of the attenuated direct 
solar irradiance and the mean intensity of the diffuse (scattered) light at the altitude, z, at solar 
zenith angle X. The attenuated solar irradiance at this height depends on the absorption by O 2 
and O 3 above this altitude, multiple scattering, and albedo. For this experiment, the 
extraterrestrial solar irradiance was not specified. However, most modelers used the spectrum 
recommended in the 1985 WMO Ozone Assessment Report, which is reported in 500 cm - 1 
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spectral intervals, the resolution used in many, but not all, of the models under consideration. 
Cross sections were in most instances taken from the most recent recommendations (JPL, 1990) 
with the exception of the O 2 Schumann-Runge band region where a number of techniques were 
employed to parameterize these values. Table K-l shows the inputs used by each of the 
participating models. 

Table K-l. Solar Irradiances and O 2 Cross Sections 


Model Solar UV O2 SRB Cross Sections P 2 Herzberg C.S. 


AER 

WMO (82) 

WMO ( 86 ) 

JPL (90) 

AFGL 

SUSIM 

Harvard-AFGL line-by-line 

Yoshino ( 88 ) 

CALJPL 

+ 

Allen & Frederick (82) 

Yoshino ( 88 ) 

CAMED 

WMO ( 86 ) 

Frederick (84) 

WMO ( 86 ) 

DUPONT 

WMO ( 86 ) 

Nicolet & Peetermans (80) 

WMO ( 86 ) 

GISS 

WMO ( 86 ) 

Fang et al. (74)* 

JPL (90) 

GSFC 

WMO ( 86 ) 

Allen & Frederick (82) 

WMO ( 86 ) 

ITALY 

WMO ( 86 ) 

Park (74) 

JPL (90) 

LLNL 

WMO ( 86 ) 

Allen & Frederick (82) 

WMO ( 86 ) 

MPI 

WMO (82) 

Allen & Frederick (82) 

WMO ( 86 ) 

NCAR 

Brasseur & Simon (81) 

Nicolet & Kennes (89) 

JPL (90) 

WASH 

WMO ( 86 ) 

Allen & Frederick (82) 

WMO ( 86 ) 

+ CALJPL solar irradiance: 




A >3300 A: WMO(1982) 

2800 < A <= 3300 A: Mentall et al. (1981) 

1172 < A <= 2800 A: Mount and Rottman (1983) 

* Includes updates for J(C> 2 ) 


The only current source of comparison to measurements is the Harvard-Air Force 
Geophysical Laboratory (AFGL) line-by-line model results, provided by Gail Anderson, which 
are based on the high-resolution spectra measured by Yoshino et al. (1983), and calibrated 
against balloon-borne high-resolution spectra (Anderson and Hall, 1986). 


Absorption by O 2 in this wavelength range (175-205 nm) is critically important to the 
photochemistry of the middle atmosphere. Because of the complex nature of this spectral region, 
atmospheric modelers have typically used analytic expressions or simplified parameterizations 
(e.g., Fang et al., 1974; Nicolet and Peetermans, 1980; Allen and Frederick, 1982). 

The following quantities were calculated as a function of z* (0-60 km) for the 
intercomparison: 


(a) column of O 2 overhead 
(c) J(0 2 ) 

(e) J(N 2 0) 

(g) J(CFC-12) 

(i) J(0 3 ->0(3P) 

(k) J(C10N0 2 ) 

(m) J(NC> 2 ) 


(b) column of O 3 overhead 
(d)J(NO) 

(f) J(CFC-ll) 

(h) J(0 3 ->0(1D)) 

(j) J(HN0 3 ) 

(1)J(n 2 o 5 ) 

(n) J(N0 3 -»N0+0 2 ) 
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RESULTS AND DISCUSSION 
(a,b) Column Densities of O2 and O3 

In principle, given that the vertical temperature structure has been specified, the O2 column 
density may be calculated in a straightforward manner. Figure K-l shows a comparison of the 
results of nine models that participated in this experiment. The AFGL results are represented by 
a solid line. With the exception of the ITALY and the NCAR models, all of the results are 
within 3% of the AFGL reference profile. Both ITALY and NCAR exhibit small variations 
about the reference with maximum differences of about 10%. 

The ozone column density is shown in Figure K-2. Somewhat larger differences are seen 
among the models, particularly near 60 km. Since the overhead abundance is dependent on the 
ozone distribution above this altitude, it is not surprising that models with differing upper 
boundary levels and specifications of the ozone abundance above the top level would exhibit 
such differences in calculated column ozone. A majority of the models are within a few percent 
of the group average. Exceptions are the ITALY and NCAR models which exhibit values which 
are systematically higher than others, by 20%-30% at 60 km decreasing towards the surface. 
Errors resulting in the interpolation of the values to the standard vertical grid may account for at 
least some of the discrepancy. 

Overall, the results of this comparison show good agreement and eliminate the column 
density calculation as a significant source of disagreement in the ensuing photodissociation rate 
intercomparison. 

(c)J( 0 2 ) 

The photolysis of molecular oxygen plays a key role in atmospheric processes. It is the 
dominant source of odd oxygen (O, Of 1 D), and O3) production in the middle atmosphere. 
However, the Schumann-Runge band region from 175 to 205 nm, which contributes much of the 
mesospheric component, is difficult to represent accurately in low resolution (e.g., 500 cm-1 
spectral interval) calculations as the absorption cross section varies by 5 orders of magnitude. 

The O2 photolysis at wavelengths less than 200 nm is depicted in Figure K-3a. A spread of 
values is seen throughout the vertical profile. In this instance, we use the high-resolution 
calculations provided by Gail Anderson at AFGL as the baseline for comparison which is 
represented as a solid line. Figure K-3b shows the ratio of the model results to the AFGL 
calculation. Those models using the family of parameterizations based on the work of Frederick 
and Hudson (1979, 1980) and Fang et al. (1974) are within 20%-30% of the AFGL results, 
particularly above 40 km. At lower altitudes, discrepancies are greater, however, the SRB region 
is of reduced importance at these levels and the accuracy of the Frederick parameterizations (i.e., 
Allen and Frederick, 1982; Frederick, 1984; WMO, 1986) is reduced at higher optical depths. 
The AER model is from 50% below to 30% above the AFGL baseline from 40-60 km. While 
this model uses the O2 cross sections recommended by Frederick in the WMO Report (1986), the 
resolution is lower than models employing 500 cm-l bins in the SRB region (approximately 1.5- 
2 nm) as it uses 48 bins of 5 nm resolution. 

Those models using other parameterizations generally show greater differences. The 
DUPONT model is systematically lower than the AFGL results by approximately 25% from 40- 
60 km. This is likely due to their use of the Nicolet and Peetermans (1980) parameterization. 
The ITALY model, using the parameterization of Park (1974), is about 30% lower than the 
AFGL results at 60 km, but is 20% higher than AFGL at 40 km, which is consistent with its 
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relatively low column O3 abundance (see Figure K-2). The NCAR model is also systematically 
lower than the AFGL calculation throughout the atmosphere. 

If the AFGL J(C>2) calculations in the SRB are taken as the most "correct" currently 
available, then the intercomparison suggests that those models using the Frederick "family" of 
parameterizations and the GISS model, which utilizes the Fang et al. (1974) SRB opacity 
parameterization, yield the closest answers to this standard. 

Recently, a polynomial coefficient representation of the Harvard-AFGL line-by-line model 
was developed by Minschwaner et al. (1992). In the future, it will be useful to employ this 
parameterization as a standard to examine the quality of more time-efficient methods for the 
calculation of O2 absorption in this region. 

When all wavelengths are included in the calculation of J(C>2), agreement between models is 
improved, as seen in Figure K-4a. Agreement with respect to the AFGL calculation is typically 
better than 15% above 30 km as seen in Figure K-4b. This improved agreement is not surprising 
as the Herzberg continuum is the dominant source of O2 absorption below 60 km and all of the 
models use recent values of these cross sections reported by the Harvard group (Yoshino et al., 
1983, 1988). Below 30 km the agreement is less satisfactory, with differences approaching 50% 
for the DUPONT, NCAR, and ITALY models. 

Calculations made with the inclusion of scattering and albedo effects show negligible impact 
above 20 km, not unexpected given that O2 absorbs UV radiation only below 242 nm. 

(d)J(NO) 

Photolysis of NO initiates the only loss of NOy in the stratosphere and therefore controls the 
roll-off of NOy mixing ratios in the middle stratosphere. 

Only the 5(0-0) at 190.9 nm and 5(0-1) bands at 182.7 nm need to be considered in the 
calculation of the photolysis of nitric oxide in the middle atmosphere. The attenuation of the 
solar irradiance in this region is due to SRB O2 absorption. Thus, results will again depend on 
the parameterization used for the SRB O2 cross sections as well as those used for the two NO 
bands. 

Figure K-5a shows the results for J(NO) for no scattering and wavelengths below 200 nm. 
This case is also representative of the other three cases owing to the location of the delta bands at 
short wavelengths. Agreement between the models is less good than that of J(02), typically 
within 30%-40% above 40 km. Both the CAMED and LLNL models are consistently lower than 
the average throughout the vertical domain, while GISS and ITALY (which do not employ the 
Allen and Frederick NO cross-section parameterization) are on the high side of the comparison. 
The ITALY model results are again well correlated with their lower than average stratospheric 
column O3 abundance. We have not yet established a reference calculation for J(NO) as we have 
for O2. 

The 30-50 km altitude range is enlarged in Figure K-5b. In this region, NO photolysis 
accounts for almost all atomic nitrogen production, which reacts with NO to destroy 
midstratospheric odd nitrogen. The spread in values is a factor of six at 30 km, decreasing at 
higher levels. 
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(e,f,g) J(N 2 0), J(CFC-ll), J(CFC-12) 

The photodissociation rates of these constituents are treated as a group because of the 
similarity of the spectral distribution of their absorption cross sections. All three constituents 
have cross sections that peak around 190 nm and fall off rapidly towards 230 nm. Because of 
this, all three should show similar behavior with respect to the SRB and Herzberg O 2 cross 
sections employed. As with J(C> 2 ) and J(NO), scattering has a negligible impact above 20 km on 
the calculated photodissociation rate. 

The photolysis of N 2 O is presented in Figure K- 6 . Above 30 km, agreement is good with all 
models within 30% of the group average. At 20 km, differences of a factor of 1.5 are evident. 

Figures K-7a and K-7b show the CFC-11 photolysis rates and their ratio to the group 
average, respectively. Better than 20% agreement is seen above 30 km, but considerably worse 
agreement in the critical lower stratosphere region. The spread at 20 km is over a factor of two 
with the AER, CAMED, GSFC, and NCAR models showing the largest deviations from the 
average. 

The CFC-12 photolysis rate and its ratio to the group average are presented in Figures K- 8 a 
and K- 8 b. The behavior is similar to that of CFC-11, with good agreement above 30 km and 
deteriorating rapidly below this level. The CAMED and NCAR models again show large 
deviations from the group average. 

(h) J(0 3 ->0(lD)) 

As the quantum yield for the photolysis of O 3 to 0( ! D) is significant only below 320 nm, 
scattering has less of a role in this rate in the stratosphere, but is still very important in the 
troposphere (i.e., Rayleigh scattering optical depth at 310 nm is about 1). Figure K-9 shows the 
photolysis rate with no scattering included. Agreement is good, particularly above 20 km with 
differences less than 10% for most models compared with the group average. The NCAR model 
shows discrepancies of up to 50% above 30 km. At 20 km factor of two differences are seen 
among all of the models. 

(i) J(0 3 -»0(3P)) 

Scattering plays an important role in the photolysis of O 3 yielding 0(3P). Figures K-10 and 
K-ll show the photolysis rate with and without the effects of scattering included reported at all 
wavelengths. With the exceptions of CAMED, NCAR, and Washington, agreement throughout 
the entire altitude range of the no-scattering case is very good with values within about 5% of the 
group average. JPL - Publication 90-1 (1990) recommends a change in the maximum quantum 
yield for the O(lD) channel from 0.90 to 0.95. The value of 0.95 was used by CAMED, WASH, 
and NCAR, thus accounting for their results lying on the low side of the group average. 

With scattering included, agreement is still very good, excepting the three models mentioned 
above, with some increased divergence in the troposphere, particularly by the GSFC model. 

(j) J(HN0 3 ) 

This molecule has an absorption spectrum below 320 nm. Once again, scattering is of 
negligible importance above 25 km. As shown in Figure K-12, agreement between the models is 
very good, particularly above 25 km. The NCAR model is consistently lower than others in the 
midstratosphere, but this discrepancy is reduced when scattering is included (not shown). This 
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may be a fortuitous result. Otherwise, agreement among the models to within 30% is seen in the 
middle and upper stratosphere. 

(k) J(C10N0 2 ) 

Figure K-13 shows the photolysis of CIONO2 with no scattering included. Agreement is at 
the 20% level above 30 km with larger differences below, excepting the NCAR and CAMED 
models, which exhibit larger differences. The CAMED model is higher than the group average 
throughout the entire vertical domain. With scattering included (not shown), enhanced rates are 
calculated below 30 km. The comparison between the models remains consistent with the no 
scattering case. 

(i)J(N 2 o 5 ) 

The N2O5 photolysis rate without scattering is shown in Figure K-14. The spread of the 
results is within 30% of the average above 40 km. The AER, GSFC, and NCAR models are on 
the low side of the other seven models presented. At 30 km, discrepancies of 60% are evident, 
increasing to a factor of two at 20 km. The comparison with scattering included is very similar 
to that of the no-scattering case. It looks as though temperature dependence is not included in the 
DUPONT, NCAR, GISS, and ITALY models, but is in the other models. 

(m) J(N0 2 ) 

With no scattering included, the agreement between models for NO2 photolysis is very good, 
to within 5% as illustrated by Figure K-15. Figure K-16 shows the result with scattering 
included. Agreement among models is within 30%. The CAMED model includes a surface 
albedo of 0.30, but does not include scattering, accounting for its low value with respect to the 
group average. 

Apparently, the inclusion of scattering leads to greater discrepancies among the models. 
Figure K-17 shows the ratio of the scattering to no scattering cases. CAMED shows a uniform 
enhancement of about 30%, which is consistent with a daytime average for a nonscattering 
atmosphere with a surface albedo (i.e., J-values increase by a factor of 1 + 2 x albedo x cos(^), 
which at noon in this example is 1.58). Other models vary from 45%-70% enhancements above 
20 km, with some models calculating nearly factor of two enhancements in the troposphere, 
while others show falloffs in this region. 

(n) J(N0 3 ->N0+0 2 ) 

Figure K-18a shows the photolysis rate of NO3 without scattering included. The models 
exhibit a considerable spread in results with differences of up to a factor of 2. Both DUPONT 
and GISS lie on the high side, while MPI lies on the low side of the other four models. These 
four models are close to the overhead sun value of 0.022 sec* 1 of Magnotta and Johnston (1980), 
which remains at the value recommended by JPL Publication 90-1 (1990). The AER model 
results are much higher than the other seven models and are not shown on the plot. 

With scattering included (Figure K-18b) the comparison among the models is substantially 
worsened. This is surprising since the wavelengths that NO3 photolyzes (around 600 nm) are 
almost insensitive to Rayleigh scattering. Only ground albedo should have a substantial impact. 
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DIURNAL AVERAGES 


A good degree of consistency is evident in the methods by which the photolysis rates are 
di urn ally averaged. Figures K-19 — K-21 show three examples of the ratio of the diurnally 
averaged rate to the noontime calculation. In each instance both Rayleigh-phase scattering and 
surface albedo are included. Figure K-19 shows the ratio for J(C>2). Only the ITALY model 
shows a significant difference from the group average, being about a factor of 1.8 below the 

raHo," 1 ^ r^ri^f 5 fe T re £ f 1116 ITALY model for a11 of the ^mal to noontime 

I !rff TH f 1 f 7 J NO) 1S P resented in R g ure K-20. In this instance, the LLNL model shows 

aiHtSc 0 p 6 V St ^ Ctl i re COn ? pared . with the other models > bein g consistently higher at all 
a titudes. Figure K-21 shows the ratio for J(N02). The curvature of the NO2 ratio in the 

troposphere is due to the Rayleigh scattering; the two models with just surface albedo and 

rTwPi 611 !? 8 at Tnl phereS (AERand CAMED ) do not show this. The increase in this ratio for 
atm : P u a “ »«* easily explained, possibly due to twilight photolysis in a spherical 
atmosphere (although this is supposedly included in some other models). See section L on 
diurnal averaging of the rates. Once again, LLNL is at variance with the group average, being on 
the low side of the other models. The tendency of the LLNL model to diverge from the grouj 

f V1 l Cn T t ^ rad P S f ° r J ( C10N °2). J(HN0 3 ), and J(N 2 0 5 ) (not shown). For all 
other constituents, the LLNL results are very close to the group average. 


CONCLUSIONS 

It is evident that a substantial improvement has been made in the consistency with which 

(JaTm 1 a S n S ^ 1 al 10I lQ r 8 a Qf Th C ^ lcala ,; ed com P ared with the previous model intercomparison 

ihn.mh'thoc!. 3 ‘Ah^ 89 • Th p 9? S .c hu ™ ann - R V n g e band region remains a potential problem, 
though those models using the family of Frederick-derived parameterizations and the Fang et al 

P“r:° t n , a T f^ eed beS } wilh , the Harvard- AFGL line-by-line calculations. At lower 
km region C ^ results are less satisfactory. Differences of 50% are seen in the 20-30 

m i depend st f on gly on Schumann-Runge band O2 absorption exhibit 

i J(N l ° ) a ^ eement 1S P° or below 40 km - The photodissociation of CFC-1 1 
and ChC-12 also show large differences in the lower stratosphere. 

CD1J Generall y S°°d agreement was seen in the photolysis rates of constituents where both the Oo 
J(HN03T° n and SCattCnng Were relatlvel y unimportant as in the case of J(C>3-*0(iD)) and 

Where scattering is important, the results are, once again, more disparate. J(C>3-^0(3p) 
discrepancies are small, though the results emphasize the importance of the O(iD) quantum yield 
parameterization. However, this difference is not important for the calculated 0(3p) since it is 
tV° ta i ? T( £» ction of O, including the quenched 0(1D), which defines its abundance (see 

fnPwr L ‘ T | 1 f e agreement . fo . r J f N ° 2 ^ 1S on ’y fair and emphasizes that the diverse approaches of 
mcludmg scattering in models lead to increased differences, as shown in Figure K-17 N?Os 

photolysis results exhibit large differences, as do those for NO3. The latter is cleariy more 
distressing in view of the recommended rate for clear sky photolysis. 
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FIGURE CAPTIONS 


Figure K-l. O 2 column abundance (molecules/cm2) as a function of altitude. Solid line 
represents the AFGL reference profile. 

Figure K-2. O3 column abundance (molecules/cm2 ) as a function of altitude. 

Figure K-3a. O2 photolysis rate (sec-1) a s a function of altitude for wavelengths less than 200 
nm and no scattering or surface albedo. Solid line represents the AFGL reference profile. 

Figure K-3b. Ratio of model O2 photolysis rates to AFGL reference profile for wavelengths 
less than 200 nm and no scattering or surface albedo. 

Figure K-4a. As in Figure K-3a except for all wavelengths. 

Figure K-4b. As in Figure K-3b except for all wavelengths. 

Figure K-5a. NO photolysis rate as a function of altitude for wavelengths less than 200 nm and 
no scattering or surface albedo. 

Figure K-5b. As in Figure K-5a except for the altitude range of 30-50 km only. 

Figure K-6. N2O photolysis rate as a function of altitude for all wavelengths with no scattering 
or surface albedo. 

Figure K-7a. CFC-11 photolysis rate as a function of altitude for all wavelengths and no 
scattering or surface albedo. 

Figure K-7b. Ratio of model CFC-1 1 photolysis rates to the group average for all wavelengths 
and no scattering or surface albedo. 

Figure K-8a. CFC-1 2 photolysis rate as a function of altitude for all wavelengths and no 
scattering or surface albedo. 

Figure K-8b. Ratio of model CFC-1 2 photolysis rates to the group average for all wavelengths 
and no scattering or surface albedo. 

Figure K-9. C>3->0(1D) photolysis rate as a function of altitude for all wavelengths and no 
scattering or surface albedo. 

Figure K-10. 03->0(3P) photolysis rate as a function of altitude for all wavelengths and no 
scattering or surface albedo. 

Figure K-ll. As in Figure K-10 except with scattering and surface albedo included. 

Figure K-12. HNO3 photolysis rate as a function of altitude for all wavelengths and with no 
scattering or surface albedo. 

Figure K-13. CIONO2 photolysis rate as a function of altitude for all wavelengths and with no 
scattering or surface albedo. 

Figure K-14. N2O5 photolysis rate as a function of altitude for all wavelengths and with no 
scattering or surface albedo. 
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Figure K-15. NO2 photolysis rate as a function of altitude for all wavelengths and with no 
scattering or surface albedo. 

Figure K-16. As in Figure K -15 except with scattering and surface albedo included. 

Figure K-17. Ratio of photolysis of NO2 for scattering to that with no scattering included. 

Figure K-18a. NO3 photolysis rate as a function of altitude for all wavelengths and with no 
scattering or surface albedo. 

Figure K-18b. As in Figure K- 18 a except with scattering and surface albedo included. 

Figure K-19. Diurnal average of O2 photolysis rate for all wavelengths with scattering and 
surface albedo included. 

Figure K-20. Diurnal average of NO photolysis rate for all wavelengths with scattering and 
surface albedo included. 

Figure K-21. Diurnal average of NO2 photolysis rate for all wavelengths with scattering and 
surface albedo included. 
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J(HN03) No Scattering, All Wavelengths 
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J(02) Diurnal Average/Scattering, All Wavelengths 
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L. Photochemistry of Radicals and Rates 

Michael J. Prather 
University of California-Irvine 


INTRODUCTION 

We examine the detailed photochemistry of the stratospheric models using the ATMOS 30N 
atmosphere to define a restricted calculation (see Table in KLM Overview section) No 
S U Nn m r S are available for such diagnostics (e.g., the 24-hour average rate of reaction ofO 
2 )> and thus we are limited to a model-model intercomparison. Nevertheless this 

budee^fo? n™ im P ortant check on the models: do we all calculate the same photochemical 
budgets for ozone given a prescribed atmosphere? The models simulated both gas-phase and 
heterogeneous sulfate-layer chemistry. We show first the gas-phase-only results and £ t^ 
relative changes induced by the prescribed heterogeneous chemistry (see Table KLM-1). 

Noontime Profiles of Radicals: O, O(lD), OH, H0 2 , NO, N0 2 , CIO. 

TJ e den fX y ° f ato T ic oxygen in da y H g ht is determined primarily by a balance between 
photolysis of O3, quenching, and recombination with 0 2 : een 


O 3 + photon 


— > 


0(3P) + 0 2 
0(lD) + 0 2 


O(lD) + (N 2 , 0 2 ) 0(3P) + (N 2 , 0 2 ) 

0(3P) + 0 2 + (N 2 , 02) -» 0 3 + (N 2 ,02) 

The lifetime of O [= 0(3P)] is much less than a minute throughout the stratosphere anH 
TCposphe^ Profiles of O and O(ID) are shown in Figures L-i anS L 2 N^ntTmeT^Uies of 
O and 0( ID) agree to within 50% with the exception of the CAMED-theta model. The CAMED 
model uses diurnal factors and calculates an effective daytime average for most of the radicals 

LLM com Pared with the noontime densities from other models. Near 40 km 

LLNL and NCAR do show some deviations from the pack, especially for OOD) The sudden 

? 5 V t 0?.Dn„ 5 ?Hn? iS inexplicab ^' M °? is the lige^ad MS 

disag^ntemVc^i^p'SSs (S s£™ "" * dircC “ y attributed to the 

(H ° x) r adicals, OH and H02, are critical to the balance of odd-nitrogen 

SSoffiM toe OHandHO H° v ™ T & mc ? hlori " e (Qy) ^mily. The system of reactions 
controlling the OH and H0 2 densities is complex and changes dramatically from the lower to 

upper stratosphere. A primary source of OH is generated by O(lD), 

0(1D) + H 2 0 -> OH + OH, 

which is augmented in the lower stratosphere by oxidation of CH4, generating 1 to 3 HOo 
radica! 5 per methane mo ecule. The individual lifetimes of OH and H0 2 are typifally less than 
100 sec, but this time scale is dominated by reactions that interchange OH for H0 2 Loss of odd 
hydrogen requires recombination of OH and H0 2 , 2 ‘ 

oh + ho 2 -> h 2 o + o 2 

or in the lower stratosphere proceeds through the NO y family, 
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OH + NO2 + (N2, O2) -» HNO3 + (N2, O2) 
OH + HNO3 -> H2O + NO3 


When heterogeneous reactions on 
stratosphere, NO2 concentrations are 
water in the sulfate particles, 


the sulfate aerosols become important in the lower 
suppressed, and the production of HNO3 and HOC 1 from 


N2O5 + (H2O, sulfate) -> 2 HNO3 + (sulfate) 

CIONO2 + (H2O, sulfate) HNO3 + HOC 1 + (sulfate). 


leads to generation of additional OH by photolysis. 


HNO3+ photon -» OH + NO2 
HOC 1 + photon — > OH + Cl 


The modeled noontime OH and H 0 2 densities (Figures L -3 and L- 4 ) agree very well above 
30 km, but diverge substantially below and are not always in proportion to the prim^ source, 
O(lD). This is particularly noticeable for the GSFC (low) and the GISS+LLNL (high)- The 
differences in the lower stratosphere apply to both OH and HO2, and thus reflect the budget for 
total HO x rather than cycling between OH and HO2. 


The NO x (NO + NO2) forms of odd nitrogen (NOv) are involved in most of the important 
chemical reactions of the NO y family. In daylight NO and N 0 2 interchange rapidly on times 
scales of 100 seconds or less, 


NO2 + photon — » NO + O 
NO + O3 — » NO2 + O2; 

however, in the lower stratosphere reactions with the Cl y family are also important in 
determining the ratio NO/NO2, 

NO + CIO -> NO2 + CI 

In daytime NO x is predicted to be typically less than 25 % of the NO y below 25 km altitude. 
Including heterogeneous chemistry reduces this to about 16 % at the minimum (about 20 km) 
under background conditions, but much less than 10% following volcanic enhancements of the 
sulfate layer. The combined NO x exceeds 50 % of total NO y above 30 km and 95 % above 40 

km. 


The modeled noontime NO and NO2 densities, Figures L -5 and L-6, are in reasonable 
agreement. The CAMED model stands out again because it calculates only "daytime averages 
rather than noontime densities. Agreement is excellent above 40 km (where all of the specified 
NOv is in the form of NO x ), but the range at 30 km and below is about a factor of 1.4 (readily 
visible in the linear scale for the NO figure). Much of this dispersion must be related to the 
relative abundance of HNO3; the predicted ratio NO/NO2 agrees to within a factor of 1.2 (see 
later discussion of Figure L- 27 ). 

Atomic chlorine (Cl) and CIO interchange rapidly through the reaction 
CI+O3 CIO + O2 
followed by 
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CIO + NO -* CI + NO 2 
CIO + O -^Cl + 02. 


Additionally, in the lower stratosphere and under volcanically perturbed conditions, models 
predict that a large fraction of CIO proceeds to Cl through the formation of intermediates such as 
CIONO2, HOC1, OCIO and BrCl. Atomic chlorine with a lifetime of less than 0.1 second is 
predicted to be a small fraction of the Cl+CIO group, and the abundance of CIO is determined 
primarily by a balance between HC1 and Cl through the following reactions: 

CI + CH4 -4 HCI + CH3 

CI + CH2O -*• HC1 + CHO 
and 

HC1 + OH -* C1 + H 2 0. 

The noontime CIO density (Figure L-7) has nearly a factor of 2 spread across the models 
throughout much of the stratosphere, but the basic shape of the profile is similar for all models. 
The CALJPL model appears to have used Cl y concentrations different from those specified in 
Table KLM-l(sce also Section M). There is some consistency in that the models with the larger 
OH (L-3) tend to have larger CIO densities. 


PHOTOCHEMICAL LOSS OF NOy 

The fall-off in NO x observed above 40 km reflects the photochemical loss of NOy. 
Photodissociation of NO in the delta bands ( see section K), 

NO + photon -» N + O 

can be followed by regeneration 

N + O2 -* NO + O 

or loss of NOy 

N + NO ->• N 2 + 0 

Figure L-8 shows that this loss rate for NO y ([N + NO] averaged over 24 hours, times 2) 
spans a factor of 4 across the models and is greatest near 40 km where the peak destruction 
occurs. This discrepancy points to long-standing problems on predicting the NO y peak mixing 
ratio near 35 km: the Harvard model (propagated here as the GISS model) always predicted very 
low maximum NOy, about 16 ppb; whereas the LLNL (not shown here) and NCAR models 
predicted more than 23 ppb. The present CAMED model predicts the least NO y loss and thus 
should have the largest NOy mixing ratios in the upper stratosphere. The GISS model has the 
largest rate of NOy destruction, but among the model results shown here it has one of the few 
independent approaches for calculating radiation in the Schumann-Runge bands (see chapter 4, 
Model Descriptions). At these altitudes, photolysis of NO occurs only in the delta-(0,0) band 
about 191 nm, and further analysis of the O2 photolysis in the S-R bands (see section K) should 
resolve these discrepancies. 
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Ozone Budget in the Upper Stratosphere: O+O3, 0+NO2, 0 +CIO, O+OH 


In the upper stratosphere the photochemical loss of O3 is dnven by reactions involving 
atomic oxygen, and we examine in Figures L -9 through L-l 2 the average rates for O+O3, O+OH, 
O+-NO2, and 0 + 00 , respectively. Focusing on 40 km altitude, we note that the spread in rates 
is much larger than the corresponding range in noontime O densities (Figure L-l). (The average 
rates from the CAMED model in the upper stratosphere fall within the range of the other models; 
such comparison provides a critical test of the CAMED diurnal averaging factors.) The model 
spread in the O+O3 rate must be due to differences in averaging over the diurnal cycle since O3 
is fixed For O+OH, the LLNL model is higher than expected from the noontime O and OH 
densities, and the CAMED model has the lowest average. The spread m the O+CIO rate is 
comparable to but larger than the range in noontime CIO densities (Figure L- 7 ), probably due to 
diurnal averaging. The NCAR model stands out here with some unusually large average rates 
above 36 km as compared with the noontime densities. 


Ozone Budget in the Lower Stratosphere: HO2+O3, HO2+NO, HOCI+hv, NO3+I1V, 
ClO+BrO, Cl202+hv 

The ozone photochemical budget in the lower stratosphere is more complex and involves 
much coupling of reactions across the H 0 X , NOv, Br y , and Cly families. Thus, the divergence o 
the model budgets is to be expected. In the lower stratosphere, the source of HO2 radicals 
involves O3. 


OH + O3 -+ HO2 + O2 

The loss of odd oxygen by this reaction (not plotted here) can be negated if followed by 


HO2 + NO -> OH + NO2, 

(shown in Figure L- 14 ) since the NO2 photolyzes to give atomic O back again. Otherwise, the 
reaction (Figure L-l 3 ) 

HO2 + O3 -» OH + O2 + O2 


or the sequence (Figure L-l 5 ) 


CI + O3 -» CIO + O2 
HO2 + CIO -+ HOCI + O2 
H 0 C 1 + photon — > OH + Cl 


results in the loss of two O3. The model calculations of these three rates agrees well with regard 
to profile shape, but still exhibits a factor of two in range. For the HO2+O3 rate, the GSFC 
model is low as expected from its H 0 X levels noted above, but the NCAR model calculates 
unusually high rates which are not consistent with its typical HO2 values. The HOC 1 photolysis 
is in excellent agreement above 25 km (except for CALJPL with different Cly) but diverges 
rapidly below. 


Photolysis of NO3 occurs rapidly in daylight ( 0.2 s*l) and about 10 % of the yield goes to O2 
(Figure L-l 6) 


NO + O3 -> NO2 + O2. 

Since most of the NO3 is formed by processes involving two O3, 
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ci + o 3 -4 CIO 

CIO + NO2 + (N 2 , 02) -4 CIONO2 + (N 2 , 02) 

CIONO2 + photon -4 Cl + NO3 

this pathway for photolysis leads to loss of two O3. 

on large ra " 8 e jr va ^ u 9 s reported for this path of NO 3 photolysis, more than a factor of 10 at 
20 km, is inexplicable in view of the previous level of agreement among the models X basic 
shape of the profiles agree, however, the divergence cannot be readily explained In the GISS 

Qopp w NO h 3 , iS f0rme ? hy P. hotol y sis of C10N02, but thVSoN02 in the GISS 

LLNL models a § ree (see Figure L- 21, and section M). Possibly, this pathway 
through NO3 is not being diagnosed in the same manner across the models. 


Halogen-driven losses of O3 occur through a coupled chlorine-bromine reaction, 

ClO + BrO -4 Br + ClOO 
■4 BrCl + O2 
-» Br + OCIO 

and the CIO dimer. 


CIO + CIO + (N2, O2) -4 CI2O2 + (N 2 , O2) 
CI2O2 + photon -4 ClOO + Cl 


The sum of the first two paths of the BrO+CIO reaction (Figure L-17) represents about 45% 
of the total reaction m the lower stratosphere and results in the loss of two O3 The dimer loss 

in , ,he An,arc, i c °r e hoie - is «* s ° 

2forB^+C10 C Fnr r,n +R S Sh0W A laT F ra ^ efor , these rates * much m ore than a factor of 
i i? T W J , For ClO+BrO the profiles have distinctly different shapes as evidenced bv the 

^FR/CA^^Wr 88 a ^ c k m (AER/CAMED/GISS vs. GSFC/LLNL/NCAR) and at 18 km 

y s ' ^ISS/LLNL/NCAR). The basic chemistry as well as the diurnal 

Hkelv t°h f /» ad,C H S i ln thC ° WCr stratos P here are more complex than in the upper stratosphere It is 
that mot 1 ^ 1 approximations of the diurnal variations in photolysis and kinetic rates is 
responsible for the growing disagreement in the lower stratosphere (see the following discussion 
tn wor?^ aVer n gmg 2‘ T hC dmrnal averaging factors used by the CAMED model do not appear 

WhCre thdr ° +N ° 2 and Cl2 ° 2+hV ratCS much 


The NO x -NO y cycling is driven by formation of HNO3 via gas-phase reaction 
OH + N 0 2 + (N 2 , 02) -4 HNO3 + (N 2 , 02), 


by heterogeneous reactions on the sulfate aerosols 

N2O5 + (H2O, sulfate) -4 2 HNO3 + (sulfate) 

CIONO2 + (H2O, sulfate) -4 HNO3 + HOC 1 + (sulfate) 

or by polar stratospheric clouds (not part of this simulation). 


HNO?f? Homing' t i on f, imposed here-background aerosol levels at midlatitudes-formation of 
HNO3 >s dominated by the gas phase reaction (Figure L- 19 ). Hydrolysis of N2O5 (Figure L- 20 ) 


contributes about 20% in the lower stratosphere. Most models predict similar values with the 
exception of CAMED (lowest N 2 O 5 from 20 to 30 km) and NCAR (unusually high N 2 O 5 at 20 
km and below). For CIONO 2 hydrolysis, the NCAR values are unusually low, and the other 
reporting models are in basic agreement (Figure L-21). These hydrolysis rates are a measure of 
the effective 24-hour average densities of N 2 O 5 and CIONO 2 . 


DIURNAL AVERAGING AND NOONTIME RATIOS 

The correct integration of radical and rates over the diurnal cycle is an essential component 
of modeling stratospheric chemistry. One part is the evaluation of photolysis rates throughout 
the daylight and twilight hours; another is the integration of the kinetic equations. Section K 
discusses photolysis rates and their integration over 24 hours. Here we examine a simple 
diagnostic of the model performance: the ratio of 24-hour- average to noontime densities of the 
key radicals O(lD), O, OH, HO 2 , and CIO. 

The density of 0( ! D) responds instantly to the ultraviolet radiation field and is created 
predominantly by the direct solar beam. Most of the models agree exceptionally well on OOD): 
In Figure L-22 the ratio increases from 0.25 in the lower stratosphere to almost 0.5 above 50 km. 
In the lower stratosphere O(lD) is formed predominantly when the sun is close to noon because 
of the large optical depths for wavelengths less than 310 nm. Above 50 km the atmosphere is 
optically thin, and O(lD) is formed at nearly a constant rate as long as the sun is up. The NCAR 
ratio is in disagreement with all the other models in this and all of these diurnal diagnostics. The 
CAMED model is a special case since they use diurnal averaging "factors." We have used their 
daytime average densities as noontime in these ratios, and thus the CAMED ratios all appear 
greater than 0.5. The interesting structure in these ratios calculated with the diurnal models 
points to the difficulties in deriving constant factors. 

The 0(3P) ratios in Figure L-23 again show similar agreement among the models except for 
NCAR. In the lower stratosphere 0(3P) is formed primarily by longer wavelength, visible light 
and hence is maintained throughout the daylight hours (i.e., ratio of about 0.5). In the upper 
stratosphere, a large fraction of 0(3p) is derived from quenching of O(iD), which is formed 
primarily about noon (i.e., ratio drops to 0.4). A probable explanation for the spread in the lower 
stratosphere is that scattered light, whose intensity varies with solar angle, is important in 
production of 0( 3 P). 

Figures L-24 and L-25 show that OH and HO 2 densities peak sharply about local noon in the 
lower stratosphere (ratios of 0 . 3 ) and become more uniform during the day at higher altitudes 
(remember, O(lD) is a primary source of HO x ). The WASH model shows some anomalies, but 
the most notable differences occur for HO 2 above 45 km where CALJPL/GISS ratios diverge 
from GSFC/WASH. It is possible that the increase in high-altitude HO 2 ratios reflects the 
treatment of the recombination of HO x after the sun sets. 

The comparison of the diurnal averaging for CIO, shown in Figure L-26, is disappointing. 
Above 45 km, CIO densities increase at twilight and at night; thus, the ratio becomes greater than 
1. The CALJPL/GISS/GSFC models agree very well, and if we accept these ratios as correct 
then the WASH model has a problem with the CIO cycle above 35 km and the NCAR values are 
systematically in error. 

The noontime ratio of NO/NO 2 shown in Figure L-27 has been a classic diagnostic of the 
chemical models. To first order it represents the ratio of the photolysis rate of NO 2 to the rate 
coefficient for the NO+O 3 rate times the O 3 density. The CAMED diagnostic refers to a ratio of 
the daytime averages and thus is expected to be less than the noontime ratios from the other 
models. Although a large portion of the disagreement can be traced to the difference in NO 2 
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photolysis rates (e.g., GISS/CALJPL have the highest rates, GSFC/WASH have the lowest; see 
section K), the remaining spread in NO/NO 2 ratios demonstrates the importance of CIO, H02, 
etc. in determining this ratio. 

Heterogeneous vs. Gas-Phase Chemistry: O, O(lD), NO, NO2, CIO, OH, HO2 

The inclusion of heterogeneous chemistry in the form of reactions of N 2 O 5 and CIONO 2 on 
the sulfate-layer aerosols (see Table KLM- 1 ) directly shifts the NO x -HN03 balance, and 
secondarily affects HO x densities and the C10-HC1 balance. For the fixed ozone and temperature 
profiles assumed here, the densities of O and O(lD) should not be affected substantially, as 
illustrated in Figures L-28 and L- 29. 

The noontime densities for both NO and NO 2 (Figures L-30 and L-31) are lowered by 
heterogeneous processing about 25% at the altitude of peak effectiveness near 18 km. The 
noontime densities for CIO (Figure L- 32), OH (L-33), and HO 2 (L- 34 ) are all enhanced by 
heterogeneous reactions: CIO by about 40%-80% at 18-20 km under these background aerosol 
levels; OH by only about 10%, but HO 2 by about 40%. The overall increase in HO x is due to the 
reduced NO x as discussed above; the larger relative increase in HO 2 is caused by the lower NO 
concentrations and the corresponding reduction in the recycling of HO 2 to OH by 

HO 2 + NO -> OH + NO 2 

The increases in noontime HO 2 are paralleled by the increases in 24-hour averaged HO 2 (Figure 
L-35). All models except for NCAR show a similar vertical pattern in these changes. The 
NCAR model is quite different and indicates very little change in noontime NO x , HO x , and CIO 
for heterogeneous chemistry; however, the diurnal average HO 2 agrees with the other results. 
(Perhaps this discrepancy is an error in diagnostics.) 

Heterogeneous vs. Gas-Phase Chemistry: Ozone Budget Rates 

The inclusion of N 2 O 5 and CIONO 2 heterogeneous processing also causes a substantial shift 
in the reactions controlling the O 3 photochemical budget in the lower stratosphere (see the 
discussion in the UNEP/WMO 1991 Ozone Assessment). A summary of the relative changes in 
the 24-hour averaged rates (see Figures L-9 through L-19 above for the gas-phase chemistry 
alone) is given in Figure L-36 through L-45. 

The average O+O 3 rate remains unchanged. As expected from the changes in the NO x , HO x , 
and CIO densities, the Cly- and HO x -catalyzed losses increase while the NO x -catalyzed losses 
become less important. The HO 2 +O 3 and NO 2 +O rates parallel the changes in HO 2 and NO 2 
respectively. The rate of O+NO 2 is reduced by 20%-40%; but HO 2 +O 3 increases by a factor of 
about 1.5, making it the dominant ozone loss mechanism in the lower stratosphere. The rate of 
HO 2 +NO is almost unchanged because of canceling effects on the HO 2 and NO densities. 
Formation of HNO 3 by gas-phase reaction (OH+NO 2 ) is suppressed slightly, but more than 
offset by heterogeneous production. 

All of the halogen catalytic cycles are enhanced by factors of 1.5 to 2 or more near 20 km. 
The NO 3 photolysis channel leading to ozone loss becomes less important in most models, but in 
the GISS model (with its enhanced CIONO 2 source of NO 3 ) it remains unchanged. Overall, the 
impact of heterogeneous chemistry on the ozone-destroying catalytic cycles in the lower 
stratosphere is consistent among the models. However, relative changes in the CIO-related 
cycles have the largest range among the models. 
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SUMMARY 


This effort in model-model comparison has shown that the stratospheric models generally 
agree on radicals and rates to within a factor of 1.5 or better. Considering the large changes in 
photochemical environment from the tropopause to the upper stratosphere, the agreement among 
models is very good. We have already identified some essential differences in model 
formulation and pointed to possible numerical errors among the suite of participating models. 
There remain many differences that cannot be ascribed a cause. Most of these discrepancies 
cannot be resolved until differences in photolysis rates are understood. 

The spread in noontime radical concentrations is worrisome, but we must resolve the 
discrepancies in J-values first. 

OH is the key radical coupling all of the chemical families that has the largest predicted range 
in the lower stratosphere; this difference should be resolved! 

The NOy loss in the upper stratosphere varies considerably from model to model; both 
photolysis and kinetic rates need to be evaluated separately. 

The diurnal cycle does not appear to be equally well represented in all models; for example 
the range in average O+O3 rate may be due to diurnal resolution as well as photolysis rates. 

The importance of NO3 photolysis in the ozone budget of the lower stratosphere is highly 
uncertain given the range in modelresults; this uncertainty carries through to the impact with 
theterogeneous chemistry included. 

Chlorine reactions leading to ozone loss show the largest spread among the models; these 
discrepancies need to be understood if we are to rely on model predictions of O3 change with 
enhanced sulfate-layer aerosols or with changing Cly levels. 

Heterogeneous chemistry has a consistent impact on the models, but is still not uniformly 
effective. 

The largest changes due to heterogeneous processing by the sulfate-layer aerosols appears in 
the 16-20 km range for which we have been and are able to make measurements of radicals 
and their cycles with the ER-2. Hopefully such data will resolve some of these 
disagreements among the models, and allow us to calibrate the impact of heterogeneous 
chemistry which is likely to be more complex than the simple parameterization used in this 
comparison. 

A review of sections K and L in a year or so might show progress and aid further model 
development, but additional ATMOS profiles (section M) are probably needed to resolve 
some differences in the lower stratosphere. 
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FIGURE CAPTIONS 

Figure L-l. O density (noon, 30N) 

Figure L-2. OpD] density (noon, 30N) 

Figure L-3. OH density (noon, 30N) 

Figure L-4. HO 2 density (noon, 30N) 

Figure L-5. NO density (noon, 30N) 

Figure L- 6 . NO 2 density (noon, 30N) 

Figure L-7. CIO density (noon, 30N) 

Figure L- 8 . 2x[N+NO] rate (24-hr avg, 30N) 

Figure L-9. [O+O 3 ] rate (24-hr avg, 30N) 

Figure L-10. [O+OH] rate (24-hr avg, 30N) 

Figure L-ll. [O+NO 2 ] rate (24-hr avg, 30N) 

Figure L-12. [O+CIO] rate (24-hr avg, 30N) 

Figure L-13. [HO 2 +O 3 ] rate (24-hr avg, 30N) 

Figure L-14. [HO 2 +NO] rate (24-hr avg, 30N) 

Figure L-15. [HOCl+hv] rate (24-hr avg, 30N) 

Figure L-16. [N 03 +hv->N 0 + 02 ] rate (24-hr avg, 30N) 
Figure L-17. [ClO+BrO] rate (24-hr avg, 30N) 

Figure L-18. [Cl 202 +hv] rate (24-hr avg, 30N) 

Figure L-19. [NO 2 +OH] rate (24-hr avg, 30N) 

Figure L-20. [N 2 O 5 +H 2 O] rate (24-hr avg, 30N) 

Figure L-21. [CIONO 2 +H 2 O] rate (24-hr avg, 30N) 
Figure L-22. 0[*D] ratio (24-hr/noon, 30N) 

Figure L-23. O ratio (24-hr/noon, 30N) 

Figure L-24. OH ratio (24-hr/noon, 30N) 

Figure L-25. HO 2 ratio (24-hr/noon, 30N) 

Figure L-26. CIO ratio (24-hr/noon, 30N) 
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Figure L-27. 
Figure L-28. 
Figure L-29. 
Figure L-30. 
Figure L-31. 
Figure L-32. 
Figure L-33. 
Figure L-34. 
Figure L-35. 
Figure L-36. 
Figure L-37. 
Figure L-38. 
Figure L-39. 
Figure L-40. 
Figure L-41. 
Figure L-42. 
Figure L-43. 
Figure L-44. 
Figure L-45. 


NO/NO2 ratio (at noon, 30N) 

O (noon) ratio (het/gas, 30N) 

OpD] (noon) ratio (het/gas, 30N) 

NO (noon) ratio (het/gas, 30N) 

NO2 (noon) ratio (het/gas, 30N) 

CIO (noon) ratio (het/gas, 30N) 

OH (noon) ratio (het/gas, 30N) 

HO2 (noon) ratio (het/gas, 30N) 

HO2 (24-hr) ratio (het/gas, 30N) 

[O+O3] (24-hr) ratio (het/gas, 30N) 

[O+NO2] (24-hr) ratio (het/gas, 30N) 

[HO2+O3] (24-hr) ratio (het/gas, 30N) 
[HO2+NO] (24-hr) ratio (het/gas, 30N) 

[O+CIO] (24-hr) ratio (het/gas, 30N) 

[HOCl+hv] (24-hr) ratio (het/gas, 30N) 
[ClO+BrO] (24-hr) ratio (het/gas, 30N) 
[Cl202+hv] (24-hr) ratio (het/gas, 30N) 
[N03+hv— >N0+02] (24-hr) ratio (het/gas, 30N) 
[OH+NO2] (24-hr) ratio (het/gas, 30N) 
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Figure L-22 (top, left) Figure L-23 (top, right) 
Figure L-24 (bottom, left) Figure L-25 (bottom, right) 
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Heterogeneous/Gas Phase Ratio Heterogeneous/Gas Phase Ratio 

Figure L-26 (top, left) Figure L-27 (top, right) 
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Heterogeneous/Gas Phase Ratio Heterogeneous/Gas Phase Ratio 

Figure L-30 (top, left) Figure L-31 (top, right) 
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Heterogeneous/Gas Phase Ratio Heterogeneous/Gas Phase Ratio 

Figure L-34 (top, left) Figure L-35 (top, right) 
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Figure L-38 (top, left) Figure L-39 (top, right) 
Figure L-40 (bottom, left) Figure L-41 (bottom, right) 






[CIO+BrO] (24-hr) ratio (het/gas, 30N) [CI202+hv] (24-hr) ratio (het/gas, 30N) 
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Heterogeneous/Gas Phase Ratio Heterogeneous/Gas Phase Ratio 

Figure L-42 (top, left) Figure L-43 (top, right) 

Figure L-44 (bottom, left) Figure L-45 (bottom, right) 
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INTRODUCTION 

The simultaneous measurement of many chemical species by the Atmospheric Trace Molecules 
Spectroscopy (ATMOS) instrument on the Space Shuttle provides a fundamental test of the 
photochemical model's partitioning of the chemical families under highly restricted conditions. 
Because the ATMOS observations define many quantities in the atmospheric profile (e.g., O3, 
H2O, CH4, NOy), this comparison tests only the photochemical components of the stratospheric 
models and is independent of many other components, such as the large-scale transport of tracers 
and predictions of ozone. 

ATMOS is a solar occultation instrument that makes nearly simultaneous measurements of NO, 
NO2, HNO3, HO2NO2, N2O5, CIONO2, and HC1. These measurements define the partitioning of 
the odd- nitrogen family (NOy) over a large range of the middle latitude stratosphere and, 
additionally, put constraints on the chlorine chemistry. We are using the ATMOS measurements of 
12 sunsets between 26N and 33N made on Spacelab 3 from 29 April to 1 May 1985. Sunset 
observations were chosen over the sunrise measurements (near 48S) because of the larger number 
and better quality of available profiles. The inclusion of sunrise data from this or subsequent 
ATMOS flights is necessary in future models and measurements comparison because the 
partitioning of radicals at sunrise provides a separate and distinct test of the photochemical models. 

The family abundances of NOy (from ATMOS), CL (from empirical correlations with NOy), 
and Bry (fixed at 15 ppt) are summarized in Table KLM-1 along with temperature, ozone, 
methane, and water vapor. The geometry for the solar zenith angles is based on a latitude of +30 
degrees with a solar declination of +15 degrees. We have chosen sulfate aerosol areas based on 
the WMO (1992) definitions for clean, non-volcanic background conditions (G. Yue, L. Poole, 
and L. Thomason, private communication). Results for both gas phase and heterogeneous sulfate 
chemistries were reported. We show only the heterogeneous results here; look at section L for 
direct comparisons of these two chemistries. 


GEOMETRY AND MODELING OF SUNSET PROFILES 

The models were asked to calculate a photochemical steady-state (i.e., repeating 24-hour 
cycles) and then to sample the species’ densities at "sunset" for the altitude range 14 to 52 km. We 
have been very careful here in defining sunset as the time of day for the tangent point along the 
occultation path from satellite to sun. This corresponds to a solar zenith angle of exactly 90 
degrees. The middle stratosphere remains illuminated by direct sunlight until the sun is more than 
4 degrees below the horizon at the surface, i.e., solar zenith angles greater than 94 degrees. 

There were many initial "bugs" in the repotting of these ATMOS profiles, primarily because 
most global models were never designed to calculate densities at a solar zenith angle of exactly 90 
degrees. The calculation requires use of a quasi-spherical atmosphere for the radiative transfer. 
The details of the rapid chemical change near sunset are of interest as a check on the chemical 
models (the reason for their inclusion in this Models and Measurements Workshop Report), but 
such details have very little impact on the overall budgets and balances, which are dominated by 
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what happens within several hours of noon. Therefore, many of the chemical parameterizations 
used in the global models were unable to compute a sunset value (e.g., the CAMED model uses 
diurnal "factors," and the NCAR model calculates 24-hour photolysis rates in deriving radicals and 
rates); others had to revise significantly their photochemical modules. In general, it was agreed 
that such comparisons are needed to test the ozone assessment models and that the photochemical 
modules in the global models should be capable of diagnosing sunset and sunrise values. 


UNCERTAINTIES OF THE ATMOS OBSERVATIONS 

Estimated uncertainties for the ATMOS/Spacelab profiles are summarized in Table M-l (see 
Rinsland, 1989; Russell, 1988; Zander, 1990, 1992). We recognize that the 12 occultation 
profiles made at sunset, ranging in latitude from 25. 6N to 32.7N over the period 29 April to 1 May 
1985, are not necessarily representative of the entire middle-latitude stratosphere. Some of the 
chemical partitioning may take days to establish, and we do not know the trajectories (i.e., 
photolytic history) of the sampled air. Only HO2NO2 and N2O5 have uncertainties large enough to 
affect the interpretation of the simulations below, although the 27% uncertainty in HNO3 
introduces a significant range to the assumed NO y profile in the lower stratosphere. Future models 
and measurements comparisons should consider that the potential systematic errors included in 
Table M-l would not affect the shape of the profile, and thus the relative profiles may be more 
tightly constrained than Table M-l implies. 


Table M-l. Uncertainties of 30N ATMOS/Spacelab 3 Volume Mixing Ratio Profiles 


Altitude p Percent Relative Uncertainty (1-sigma) 


(km) (mb) 


51.5 

0.7 

47.5 

1.2 

43.5 

1.9 

39.5 

3.2 

35.5 

5.5 

31.5 

9.7 

27.5 

17 

23.5 

32 

19.5 

60 


NONO 2 


20 

25 

20 

25 

20 

15 

20 

15 

20 

15 

20 

15 

20 

15 

20 

15 

20 

15 


HNO 3 


HO2NO2 N7O5 CIONO2 HC1 
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10 

27 

70 
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10 

27 

40 

32 

12 

10 

27 

40 


15 

10 

27 
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NO2 and NO [Figures M-l and M-2] 

The agreement of ATMOS measurements and the models for NO and NO2 is surprisingly 
good. In the upper stratosphere, these species comprise the majority of the NO y , and thus the 
specification of NO y forces some of this agreement. However, the fall-off of NO2 above 35 km is 
a true test of the models, and it is excellently reproduced. In the lower-middle stratosphere, given 
the number of factors affecting NO x /NO y partitioning (e.g., OH, heterogeneous chemistry, three- 
body rates), the model agreement with ATMOS for both NO and NO2 is a very successful test. 
Differences over the critical region 16-26 km are less than 30%. There is a systematic tendency for 
the models to underpredict the NO2/NO ratio by about 20% throughout the profile, but one must 
consider the difficult geometries in calculating NO2 photolysis through sunset. 

Grade: A- 
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HN 0 3 [Figure M-3] 


The models agree with the ATMOS profile for HNO 3 to better than 20% over most of the 
stratosphere. A most important region is the fall-off above 30 km where the HNO 3 abundance 
reflects a balance between the OH + NO 2 reaction and photolysis. The comparison here, where 
HNO 3 is a small fraction of the NO y family, is excellent. 

Grade: A 

HO2NO2 [Figure M-4] 

The majority of the models (CALJPL, GISS, LLNL) agree with the ATMOS profile over the 
limited range (20-34 km) to better than 30%. The GSFC model is generally larger by as much as a 
factor of about 2. The shape of the profile with a peak mixing ratio near 28 km is accurately 
predicted by all models. 

Average Grade: B+ 

N2O5 [Figure M-5] 

The difference among models, a factor of 2, is larger than for the previous NO y species. The 
ATMOS measurements are systematically smaller than all of the model predictions by at least a 
factor of 3. At sunset, N 2 O 5 has been photolyzed all day and is predicted (and observed) to be a 
small fraction of the total NO y species; near sunrise, however, there is an accumulation of N 2 O 5 
resulting from its formation throughout the night. (In the GISS model the ratio of sunrise to sunset 
is about a factor of 6 at 30 km and smaller below 30 km.) An earlier model analysis of the 
ATMOS N 2 O 5 measurements (Rinsland et al., 1989) showed that the sunrise data (taken at 48S, 
not part of the current comparison) are in reasonable agreement with model predictions. One 
possible cause of the discrepancy shown here might be that photolytic loss of N 2 O 5 during the day 
is underestimated by the models (see discussion in Rinsland et al., 1989); however, the 48S 
profiles represent quite different atmospheric conditions in terms of temperature and photolysis 
rates, and the combination of 48S at sunrise followed by 30N at sunset cannot be viewed simply 
as a measure of the diurnal cycle of N 2 O 5 . 

Nevertheless, the overestimate of N 2 O 5 is worrisome, because our heterogeneous chemical 
removal of NO x to HNO 3 , proportional to N 2 O 5 , would also be in error. (However, the 
NOx/HN 0 3 , ratios noted above look good.) It is important to re-evaluate this problem, focusing 
on altitudes below 25 km in light of the importance of the N 20 s-sulfate reaction. We recommend 
extending this comparison to the 48S sunrise data and to the recent 1992 ATMOS measurements. 

Grade: C 

CION 0 2 [Figure M- 6 ] 

The model predictions are in good agreement with each other where the peak levels of CIONO 2 
occur between 26 and 34 km. However, they are uniformly greater than the ATMOS profiles by 
about a factor of 1.5. An exception is the CALJPL model, which shows excellent agreement 
between 20 and 30 km; however, the CALJPL model chose a Cl y profile, different from that 
recommended in Table KLM-1, based on their previous analysis of the ATMOS profiles. (The 
CALJPL choice appears optimized to CIONO 2 rather than HC1; see next section.) Possibly, some 
of the problems with overpredicting N 2 O 5 at sunset apply to CIONO 2 , which is also photolyzed at 
long ultraviolet wavelengths. Overall, the shape of the profile is excellent. 
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Grade: B+ 


HCI [Figure M-7] 

The predicted HCI profile shows a lot of structure in the lower stratosphere, with a minimum 
mixing ratio in some cases near 28 km at the CIONO2 maximum. The range in model results for 
HCI between 20 and 40 km can in general be attributed to the different predictions for CIONO2 
below 30 km and that for CIO above. The ATMOS profile does indeed show an inflection point 
near 24 km, but its shape is different from the models. One source of this discrepancy may be the 
prescribed Cl y profile (based loosely on aircraft and balloon data correlations of halocarbons with 
N2O and thence NO y ), which may be too large in the lower stratosphere. (The CALJPL profile 
has less Cl y in the lower stratosphere.) There still remain some obvious problems with the Cl y 
partitioning, because if we look at the ratio CIONO2/HCI all of the models (including CAUPL) 
would be systematically higher than ATMOS. 

Grade: B- 

CIO and HOCI [Figures M-8 and M-9] 

We include the modeled sunset profiles of CIO and HOCI for completeness, although ATMOS 
did not measure them. Although all of the reporting models produce a similar profile for both CIO 
and HOCI, there are some large differences. CALJPL has, of course, lower CIO mixing ratios 
between 20 and 30 km because of their lower Cl y profile; however, that does not explain the much 
larger range of results for HOCI. The GSFC model diverges above 40 km, predicting substantially 
more CIO and HOCI than the other models. However, the models (AER, CALJPL, GISS, GSFC, 
LLNL) agree reasonably for HOCI and CIO (and hence Cl y partitioning) in the lower stratosphere 
near 20 km. In future models and measurements comparisons it may be possible to use the 
simultaneous measurements of CIO and HCI from the ER-2 to augment the ATMOS test in the 
lower stratosphere. 

BrO, BrON0 2 , and HBr [Figures M-10, M-ll, and M-12] 

The predicted partitioning of key species of the bromine family are also included here for 
completeness of documenting the models. The sunset values of BiO range from 2 to 4 ppt near 20 
km where the bromine catalyzed ozone loss is important. Throughout the lower stratosphere 
B1ONO2 is expected to be about 10 ppt and is the dominant form of bromine. (The LLNL model 
has an unusual B1ONO2 profile below 25 km and must have the Br y in some other form.) HBr is 
expected to be small, less than 1 ppt, in all reporting models; but the GSFC model here is aberrant. 
It is not clear whether we will ever have any pair of Br y species measured that would test this 
partitioning. 


SUMMARY 

In summary, the ATMOS comparison was surprisingly good once the "sunset" reporting 
difficulties were overcome. The N2O5 discrepancy is worrisome, but the NO x /HN03 partitioning 
in the lower stratosphere looks good. More effort needs to be put in the Cl y analysis, perhaps 
resetting the profile, before we can examine the HCI profiles in detail. 

Despite the large range (factor of 1.5 to 2) in model profiles relative to the ATMOS 
measurements, this comparison is quite encouraging. For example, most of the differences 
correspond to less than a 2-km shift in the profiles. If we are incorrectly modeling any key 
reactions that control such partitioning then it must have a profile similar to the currently modeled 
reactions. Therefore, matching the observed ATMOS profiles over two orders of magnitude 
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throughout the stratosphere demonstrates that the basic mechanisms partitioning the chemical 

families are probably well represented in the models. 
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FIGURE CAPTIONS 


Figure M-l. NO 2 (het, sunset, May 30N). 
Figure M-2. NO (het, sunset, May 30N). 

Figure M-3. HNO 3 (het, sunset. May 30N). 
Figure M-4. HO 2 NO 2 (het, sunset. May 30N). 
Figure M-5. N 2 O 5 (het, sunset, May 30N). 
Figure M-6. CIONO 2 (het, sunset. May 30N). 
Figure M-7. HC1 (het, sunset, May 30N). 
Figure M-8. CIO (het, sunset, May 30N). 
Figure M-9. HOC1 (het, sunset, May 30N). 
Figure M-10. BrO (het, sunset, May 30N). 
Figure M-ll. B 1 ONO 2 (het, sunset, May 30N). 
Figure M-12. HBr (het, sunset. May 30N). 
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INTRODUCTION 

riiv J rtle ,ocal 1 tendency of the continuity equation can be separated into terms representing transport 
p ocesses and photochemical production/loss. In order to understand the modeled response ofthe 
trace gases to changes in the photochemical environment, it is desirable to compare ^magnitudes 

timescale 61 ™ 5 ' ^ expenment re P resents a fast effort at examining the balance on a seasonal 

. . Jt 1S cl j ? ar that , the balance depends on the distribution of the trace gas being considered. For 
instance, the steady-state mixing ratio for an inert tracer is constant in space and time, and the 
transport fluxes due to advection and eddy diffusion vanish at every time-step. In this experiment 
ave examined O3 and N2O. Nitrous oxide has also been chosen to illustrate a typical upward- 

seasmvaf aim 1S t0 der * ve a methodology to illustrate which terms determine the 

sonal tendency of the trace gas. A next step would be to perform appropriate spatial and 
temporal averages to understand the year-to-year and steady-state balance among the terms 


THE EQUATION 

The continuity equation can be written in the form 


3f . 

= +S a dv + Sdiff + P - L, 

t ? nden e y tenn { u r the ™ ixing rati0 f; S adv is the tendency due to transport by the 
circulation, Sdiff is the tendency due to eddy transport parameterized by the eddy 

removarrat?^ 16111 ’ ? 1S ^ 003 photochemical production rate; and L is the local photochemical 


The terms S a dv and Sdiff can be separated into the horizontal and vertical components, 
exact forms tor S a dv and Sdiff depend on the formulation and coordinates being used, 
example, m pressure coordinates, the term S a dv can be written as 


The 

For 


c 9 (vf) 3 (wf) 

^ adv ~ 5z~ 


or 


df df 
V ~' w di 


dy 


th * ? add , up to the same sum in each case, the separation into the horizontal and 

vemcnl components is different depending on whether the differentials are calculated as finite 
differences between adjacent grid points or directly as fluxes across boundaries of the grid boxes. 

For Sdiff, 


Sdiff- -^( Kyy - ^(Kyz^) 

<Lir K 

- & Kz V ’ & Kzz P * 

In displaying the results, we have used the following labels : 
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H a dv = ■ [flux form] or 
V adv = *^^ [flux form] or 


df 

[gradient form] 


df 


- [gradient form] 


Sadv = Hadv + ^adv 
Hdiff= - j^ K yy||) - 

Vdiff= - ^(Kzy^) - ^ Kzz ^ } 


Sdiff = Hdiff + Vdiff 


TRAN = Sadv + Sdiff 
TEND = ^ = TRAN + P - L . 


BRIEF DESCRIPTIONS OF THE MODELS 

Five models contributed results for this experiment. Table N-l summtmzes some ^of the key 
differences among the models. For more detailed descnptions of the model formulation, see 

chapter 4. 


Table N-l. Model Descriptions 


Model 

Coordinate and 

Numerical 

Scheme 

Circulation 

Eddy Diffusion in the Stratosphere 

AER 

log-pressure; 

Smolarkiewicz 

diabatic circulation; 
prescribed 

Kyy = 3x109 -lxlOWcm 2 sec- 1 , 

Kzz= 1x103 cm2 sec- 1 , 

Kyz by projection of isentropic surfaces on 
to pressure surface 

CAMED- 

6 

isentropic; 

Adam- 

Bashforth 

interactive 

Kyy, calculated, different for N 2 O and O 3 
Kzz = 3x103 cm 2 sec 1 , 
no Kvz 

GSFC 

log-pressure; 
Prather scheme 

residual mean from 
prescribed diabatic 
heating 

Kyy, calculated, 

Kzz = 2x103 cm 2 sec* 1 , 

Kyz by projection of isentropic surfaces on 
to pressure surface 
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ITALY 

T T XTT 

log-pressure; 

centered- 
differencing 
with Lorentz N- 
cycle time 

integration 

diabatic; 

from 3-D spectral 
model 

Kyy calculated from 3-D spectral model, 
typical values 1-5x108 cm2 sec-1, in the 
tropics, 3-8x109 cm2 sec-1 extra-tropics 
Kzz < 5x102 cm2 sec-l tropics, 103 cm2 
sec-l in the extra-tropics 
Kyz by projection of isentropic surfaces on 
1 to pressure surface 

f 1-/NL/ 

1 log-pressure; 
Smolarkiewicz 

residual mean from 
prescribed diabatic 
heating 

Kyy, 2x109 cm2 sec-l, 

Kzz = 1x103 cm2 sec-l, 

Kyz by projection of isentropic surfaces on 
to pressure surface 


“iss sasasa. ns s - •" - 

DATA COLLECTION 

lorn'S Sd£ SeSr by the da,a base ' ^ foll °-"S 'WO points 

(1> S ! .' h ! l, r fc differem fonnulations in the models, some of the quantities requested mav 
*f. fo™ “sed m the actual time integration of the model equations [n those ca^ 
the quantities are calculated off-line as diagnostic variables Thn« • * y 

quantities do no, necessanly explatn diffemSmLTltrul.s ' * ,heSe 

(2) The data requested represent a snapshot on a particular day. While it has been verified in 
some of the models that the snapshot is representative of die monthly mean this has nor 

the vai° ne f ° r a thC models ‘ Des P ite this, the discussion in this section will assume that 
the values are representative of the monthly mean. 


RESULTS FOR N 2 0 

KpsSs sBS: 
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Table N-2. The Rate of Change Thai Would Affect the Given Mixing Ratio by 20% Over a 90 



The following points are noted : 

m There is no in situ production of N 2 0 in the atmosphere. It is transported from the 
1 troposphere to the stratosphere where it is removed by P ho ^ h emical mcim N 2 C )s 
vertically stratified tracer, i.e., the mixing ratio decreaseswithheight. All modells ishow ^a 
downward sloping of the surface of constant mixing ratio from the tropics towa d 

pole. 

stratosphere. This would imply significant differences m the seasonal behavio 2 
those regions if the values given are interpreted as seasonal changes. 

"T-Up, c term anc i the Sdiff term are opposite in sign. S a dv is positive in the tropics and 
negative at high latitudes, consistent with the upward motions and downward motions in 
the respective regions. The behavior of S d iff (negative in the tropics and positive at h gh 
latitudes') is consistent with the down-gradient transport expected from the eddy-diffusion 
StoT Magnitude of S a dv is larger than .ha, of S ditr in the troptcs. resulttng tn a 

positive TRAN term. 

(A) In a ii t he models that use the flux form, the V a dv term carries the N 2 0 from the lower 
troDical stratosphere to the mid stratosphere to maintain the concentration against 
nhotochemical removal. Both H adv and H d iff help cany material from the tropics to he 
extra tronics The downward motion in the extra tropics also shows up as a negative 
V H The behavior of H adv in the CAMED-0 model is different because the 
tai. ^ted aTTSi Uttfl* in the same direction, H adv is everywhere 

positive. 

a general comment is that the fields in the GSFC and LLNL results show more vertically 
(5) layeredltructure, while the AER, CAMED-0, and ITALY model results are smoother. 
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The SAMS monthly mean data have limited utility in providing values for the observed 
tendency. Jones and Pyle (1984) estimated uncertainties of -15% or 30 ppbv for the derived 

!hp n t?nH mCan mix jng ratio in the lower stratosphere. This corresponds to an uncertainty in 
the tendency of about 10-14 volume mixing ratio sec-1. Both the 1979 and 1980 data show 

^40 km L^ToT ^ X 1015 V , 0lU T mking rati ° «■' ±10 * iO' 15 volume 

fhp 4 ? ^ and 30 k T’ ^P ectlvel y- However, the spatial behaviors of the tendencies from each of 
the 2 years are quite different . More data are needed to obtain climatological mean values with 
smaller uncertainty for comparison with model results. 


RESULTS FOR OZONE 


The results for 0 3 are shown in Figures N-6 to N-10. The (a) panels of Figure N-6 to N-10 
show the calculated mixing ratio TCND, TRAN, and P-L. The mixing ratio is given in ppm v All 

th!t n of change m 10-15 volume mixing ratio sec-l. The following table relates the rate 
that will give rise to significant seasonal changes for a given mixing ratio. 


Day Period ^ ° f Change That Would Affect the Given Mixing Ratio by 20% Over a 90- 


Local Mixing Ratio for O3 (ppmv) 

kate (in llH^ sec-l ) That Will Begin to Affect 
Seasonal Behavior 

10 

"260 

8 

"2l0 ” 

6 

150 

4 

loo 

2 

l0 

1 ~ H 

“26 

0.5 

13 

0.2 

1 ~ 

0.1 

13 


imrirtln? ° f ° 3 where there ar e both in situ photochemical production and removal, it is 
portant that one compare the transport terms with the P term and the L term separately instead of 
just looking at the net term (P-L) (see Ko et al., 1989). We added a set of figures for the ozone 

Md1te 1 q n u ^tiv E^fll2d^s 8 “ reS N ' 6 thTOU6h N ' 10 ’ ' VWCh Sh ° W ,he individual P a " d L. 

PO - a MP) - abs(L) 
abs(P) + abs(L) 

where abs(.) is the absolute value. Noting that when P is equal to L in the photochemical 
equilibrium region, EQ = 0 in that region. The condition EQ > 0 marks theregion where 
production is balanced by transport away from the region, while EQ < 0 marks the region where 
local removal is balanced by transport into the region. * 

The following comments can be made concerning the O3 results: 

• Above 30 km in the tropics, both the P and L terms are at least two orders of magnitude 
larger than any of the transport terms. Thus ozone is photochemically controlled in that 
region as is evident from the top panel of (c) in each model. 
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. In the lower stratosphere, production in the tropics is balanced by transport away from the 
region °For results presented in flux form, both V adv and H adv play a role m exporting the 
ozone from the lower tropical stratosphere. In most of the models, the magnitude of S adv is 
teger ton iS It is an oZ of magnitude larger in the trop.cs and a factor of 2 

larger in the extra tropics. 

Tr» further examine the seasonal behavior of ozone, we present a number of figures 
representing the local lifetimes defined by each process. We define the lifetimes corresponding o 

each tendency term as follows : 


f 

t PROD - a bs(PROD) 

f 

t TEND- abs(TEND) 

f 

T S a dv- a bs(S adv ) 


f 

t LOSS - abs(LOSS) 
f 

t tran- abs(TRAN) 
T S d iff- a bs(S d iff) 


are shown in Figures N-l lto N-15. 


We note that : 

Thp t in the tropics around 20 km is of the order of 1000 days for all seasons, 

c(?nsistent N with the lack of observed seasonal variation in that region. The differences among 
the models are probably insignificant. 

. In the tropics between 16 to 24 km, photochemical production is balanced by TRAN. 
Ab?ve tomodel tropopause, S a(iv is larger too Sup However, to behav.or of s adv and 
Sdiff around the tropopause is very different among the models. 

• In the tropics between 16 to 24 km, all model results show that the lifetime due to 
photochemical removal (L) is about 1000 days, much longer than the lifetime due to P or 
TRAN h This would suggest that ozone in the tropics would not be very sensitive to chemical 
nfrtu^tions There M some differences in the P and TRAN. The model results separate 
P „ tu p &ER GSFC and ITALY results indicate that the lifetimes due to P 

SflXSti! ta iffi. t to cases of LLNLand CAMED-9, to lifetimes are 

shorter than 100 days. 

. The seasonal signal is much stronger in the mid- and high-latitudes. For March, the net 
mndenc v°Ufetime at z* = 20 km is 1000 days and longer, with advection balancing eddy 
HiffnsS and local loss. In June most models show that all the terms contributing to the 
TRAN are negative, leaving a negative tendency with a lifetime of about 200 to 500 days. 
The exception is GSFC, which shows a positive tendency due to advection, resisting l 
positive TRAN term. In December there is a positive net tendency for all models with 
lifetimes of 200 to 500 days where both the advection and diffusion terms are positive. 
A^aim ^he exception i s in the GSFC case in which the diffusion term is negative. 

A -r of 200 days implies a change of 15% in one month. A more typical value of 400 
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FIGURE CAPTIONS 


Figure N-l Model calculated results for N20 for March 15 from AER model. 

The contour intervals are: ±10, ±5, ±2,±i,2W.o, xu.i, ™ » 
of the labels, (b) Six panels — individual tendency terms for N 2 • 

Figure N-2. Same as Figure N-l except the results are from the CAMED-0 model. 

Figure N-3. Same as Figure N-l except the results are from the GSFC model. 

Figure N-4. Same as Figure N-l except the results are from the Italy model. 

Figure N-5. Same as Figure N-l except the results are from the LLNL model. 

the quantity EQ (see text). The quantity EQ is dimensionless. 

Figure N-7. Same as Figure N-6 except the results are from the CAMED-0 model. 

Figure N-8. Same as Figure N-6 except the results are from the GSFC model. 

Figure N-9. Same as Figure N-6 except the results are from the ITALY model. 

Figure N-10. Same as Figure N-6 except the results are from the LLNL model. 

Figure N-ll. Model calculated lifetimes for specific mechanism from the AER model for 

mLh. June and December oftSa The 

with a scud Hue. If .he iendency is 

negative, t is plotted with a dashed line. 


Figure N-12. Same as Figure N-ll except the results are from the CAMED-0 model. 
Figure N-13. Same as Figure N-ll except the results are from the GSFC model. 
Figure N-14. Same as Figure N-l 1 except the results are from the ITALY model. 
Figure N-15. Same as Figure N-ll except the results are from the LLNL model. 
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Figure N-9b 
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INTRODUCTION 

These simple experiments allow us to calibrate the transport of the different models, isolating 
the effect of transport on lifetimes and tracer distributions. We define two tracers, XI and X2, 
which have their mixing ratio set to 1 ppbv everywhere below 500 mbar. The loss frequency, 
L(p), for XI and X2 is only a function of pressure p, is independent of latitude and season, and is 
defined as the following: 


L(p) for XI 

L(p) for X2 

Pressure range 

(sec- 1 ) 

(sec- 1 ) 


0.0 

0.0 

p > 100 mbar 

3.0E-6/p2 

1.5E-6/p2 

1 < p < 100 mbar 

3.0E-6 

1.5E-6 

p < 1 mbar 


Modelers were asked to run the experiment to a steady-state distribution and report the steady- 
state lifetimes and the latitude-z* grid of monthly mean mixing ratios for XI and X2. 

The loss frequencies defined for XI and X2 were picked to give lifetimes for these gases that 
were around a century or more. These lifetimes are similar to those expected for source gases N2O 
and several of the chlorofluorocarbons. The results of these idealized model studies are, therefore, 
relevant to sections E (Large-scale Structures in N2O and CH4) and H (Correlation of Long-lived 
Species in Simultaneous Observations). Since the distributions of XI and X2 are dependent on the 
transport of the individual models, we rely on some information sent for section A (Temperatures 
and Net Radiative Heating) in the analysis of these results, namely, the vertical velocity. 


MODEL RESULTS AND ANALYSIS 

Tracer X2, with the smaller loss frequency, is transported to higher altitudes and thus fills 
more of the 0 to 60 km domain of the UADP grid. We have, therefore, plotted the March, June, 
September, and December latitude-z* grids of X2 for each model in Figures 0-1 to 0-10. The 
distribution of tracer X2 is affected by the total transport (both advection and diffusion); however, 
there is a strong correlation with vertical velocity, w, which is also plotted in Figure 0-1. The 
influence of vertical advection on X2 is clearly evident for all models. Strong regions of upward 
motion apparently control the height and latitude of maximum vertical extent of X2 in all models. 
The mixing ratio contours are clearly responding to the circulation changes over the course of a 
year (see all plots in Figures 0-1 to 0-10). The differences in tracer distribution and vertical wind 
fields are also apparent among the 10 models shown here. 

A lifetime of XI (or X2), Lxi (or Lx2). is defined by the ratio of the total burden over the total 
loss rate, all averaged over a year. If we maintain a steady flux of XI (or X2) such that the 
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concentrations exactly repeat an annual cycle, then we derive a steady-state lifetime. If the flux is 
stopped, then the annual cycle of concentrations will decay from year to year with a 1/e folding 
time that is very close, but not identical, to the steady-state lifetime. (The 1/e folding time in this 
transient decay period will be slightly faster than the steady state because the stratosphere, which 
lags behind the troposphere, will have a slight excess of XI (or X2) compared with steady state 
and hence slightly greater loss.) The lifetimes Lxi and L\2 are given in Table 0-1. 


Table 0-1. Lifetimes of XI and X2 and Their Ratio for the Individual Models 


Model Lxi(years) L\ 2 (years) Lxi/Lx2 


AER 

117 

146 

0.801 

CALJPL 

101 

124 

0.815 

CAMED 

114 

145 

0.786 

DUPONT 

122 

152 

0.803 

GSFC 

98 

134 

0.731 

ITALY 

110 

142 

0.775 

LLNL 

105 

132 

0.795 

NCAR 

101 

134 

0.754 

NOCAR 

92 

118 

0.780 

WASH 

127 

161 

0.789 


The range of Lxi is from 92-127 years and of Lx2 is from 118-161 years. The ratio of the 
largest to the smallest lifetime for XI is 1.38 and for X2 is 1.36. The range of Lxi/Lx 2 is from 
0.731 to 0.815 and the ratio of the largest to the smallest value of Lxi/Lx 2 is 1.11. Generally, a 
short lifetime for XI implies a short lifetime for X2 when compared with other models. 

The models can be grouped according to XI lifetimes as follows: 

Short lifetimes - CALJPL, GSFC, NCAR, and NOCAR 

Medium lifetimes - CAMED, ITALY, and LLNL 
Long lifetimes - AER, DUPONT, and WASH. 

A comparison of our groupings of models and their respective lifetimes to those for carbon- 14 
(in Table I-2e of section I) show some similarities and some differences. There exists agreement in 
groupings of the models between this classification and those in section I for GSFC, NCAR, 
CAMED, LLNL, DUPONT, and WASH. Models AER, CALJPL, and ITALY are grouped 
differently in the two sections. Since carbon- 14 is input at high latitudes in the lower stratosphere 
and XI is input at the ground, a straightforward total consensus in lifetime groupings between this 
section and section I is not expected. Generally, however, a model that computes a shorter lifetime 
for XI will tend to predict a shorter residence time for carbon- 14 and vice versa. 

The lifetimes are mainly determined by the speed at which the vertical advection deposits the 
tracer in its primary loss region. For example, a large upward motion in the equatorial region 
should lead to more tracer XI in the middle stratosphere, which will lead to more XI loss and a 
shorter lifetime. We show plots of tracer XI as a function of month for several latitudes and 
altitudes in Figure 0-11. Models with short Lxi's show larger abundances of XI at 20 and 30 km 
near the equator and 35N (see NOCAR and GSFC values in Figures 0-1 la-d) than those models 
with long Lxi's (see WASH and DUPONT values in Figures 0-1 la-d). At higher latitudes, it is 
not clear which model should show the larger XI abundance. A large upward vertical motion in 
the tropics must be compensated by large downward vertical motions at other latitudes due to the 
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balance of mass motion required across a pressure level. Thus the X 1 abundance at high latitudes 
for models with short lifetimes may be depressed relative to the other models. The reverse may 
also be true: the XI abundance at high latitudes for models with long lifetimes may be enhanced, 
relative to the other models. Such is the case of WASH at 70N and 30 km (Figure 0-1 If). 

The seasonal change in the transport is also apparent in Figure 0-11 for latitudes 35N and 
70N. In Figure 0-1 Id most models show a seasonal dependence with maximum mixing ratios of 
tracer XI predicted in the late summer and early fall time period; presumably occurring after a few 
months of upward motion. The only model that shows a different pattern from this behavior is 
LLNL (with maximum XI abundances in late winter). Although all models at 70N (Figures 0-1 le 
and 0-1 If) show generally the same pattern with maximum XI abundances in the summer to fall 
months, there are model differences in predicting the peak X 1 abundance (e.g., ITALY shows a 
peak in June or July, while LLNL shows a peak in November at 30 km). 

Plots of X2 are shown in Figure 0-12 for all models. We show the same latitudes and 
altitudes as shown in Figure 0-11, but also extend to 40 km because of the smaller loss frequency 
of X2. We expect that a large upward motion in the equatorial region should lead to more tracer 
X2 in the middle stratosphere, which will lead to more X2 loss and a shorter lifetime. This 
statement is more difficult to prove with X2 because the smaller loss frequency of X2 allows a 
larger abundance of X2 into more regions of the model regime; therefore, the X2 abundance is 
controlled by advection over a larger region. NOCAR (with the shortest Lx2) shows the largest 
X2 abundance at 20 km, approximately the third largest X2 abundance at 30 km, and about the 
second largest X2 abundance at 40 km (Figures 0-12.a-c). Clearly the X2 abundance at 35N also 
makes a difference in the calculation of Lx2 (notice that NOCAR is at or near the top in X2 
abundance for 20, 30, and 40 km in Figures 0-12d-f). WASH (with the longest X2 lifetime) 
fairly consistently shows abundance of X2 at or near the bottom of all the models at both the 
equator and 35N (see Figures 0-12a-f). 

A model with strong upward motion in one portion of the stratosphere does not necessarily 
have strong upward motion throughout the entire stratosphere. GSFC, with the second shortest 
Lxi, has the fourth shortest Lx2- Both CALJPL and LLNL have longer lifetimes than GSFC for 
XI, but have shorter lifetimes for X2. GSFC shows the second largest X2 abundance at the 
equator and 30 km (Figure 0-12b), but at the equator and 40 km (Figure 0-12c) GSFC has about 
the smallest X2 abundance. Clearly the upward motion is slowing down significantly (when 
compared to other models) in GSFC between 30 and 40 km. 

Much of the seasonal behavior in the X2 abundance from a particular model is similar to that 
pointed out for XL It is interesting to note that at 35N and 40 km (Figure 0-12f) the LLNL model 
shows similar behavior to other models with a peak in the late summertime, while at 35N and 30 
km (Figure 0-1 2e) the LLNL model shows a very different behavior compared with the other 
models. The WASH model shows a very substantial equatorial variation at 20 km in both XI and 
X2 that is not predicted in other models (Figures 0-1 la and 0-1 2a). 

The altitude behavior of XI and X2 is given at five latitudes (70S, 35S, O, 35N, and 70N) for 
March in Figures 0-13 and 0-14. This allows for inspection of the very different altitude structure 
of the tracers from the various model simulations. No one model shows the largest abundance of 
XI or X2 throughout the entire altitude range for any of the latitudes. The models (as noted 
earlier) are each different and unique. The movement of XI and X2 in each model cannot be 
adequately described in this brief report, but we do mention some of the more interesting model 
behavior below. 

For example, the ITALY model has nearly the smallest abundance of XI at 35S (Figure O- 
13b), but the largest abundance of XI at the equator (up to about 30 km; Figure 0-1 3c). The 
WASH model predicts the smallest abundance of XI up to about 23 km, whereas above 28 km the 
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WASH model predicts the largest abundance of XI for a few kilometers (Figure O- 13a). Clearly 
substantial transport of XI (probably from lower latitudes) is taking place in the WASH model that 
causes this behavior at 70S. Some models show substantial upward transport in certain regions, 
such as GSFC between 20 and 40 km at 35S (Figure 0-1 3b) and LLNL between 25 and 35 km at 
35N (Figure 0-1 3d). 

Long-lived species should exhibit approximately linear interrelationships in regions of the 
atmosphere where the lifetimes are long compared with mixing (Plumb and Ko, 1992). We 
demonstrate this statement with plots of XI versus X2 in Figure 0-15. The model results are 
indicated by circles, squares, and triangles for different latitude ranges. The solid line is the 
theoretical slope (i.e., the ratio of upward fluxes) and indicates the ratio of Lxi to Lx2 for each 
model. All the models show a very similar shape with a ramp up to the solid line and approximate 
agreement with the solid line at values of XI above 0.3 ppbv and X2 above 0.4 ppbv. Some 
models have more scatter (GSFC, ITALY, and NOCAR), while one model (CALJPL) has very 
little scatter in the plotted points. 


TRACER W 

Another tracer, W, was simulated by the AER, GSFC, and WASH models for a separate study 
connected with the NASA "Report on Concentrations, Lifetimes, and Trends of 
Chlorofluorocarbon s, Halons, and Related Molecules in the Atmosphere." Tracer W has a loss 
frequency equal to that of XI but has a temporally changing input flux. The W flux is increased 
from zero over a course of 10 years, stays constant for a decade, and then is decreased to zero by 
year 30 (as shown in Table 0-2). All model simulations were for 40 years. Tracer W was chosen 
to be similar to a chlorofluorocarbon, with fluxes increasing in the 1960s and 1970s, leveling off 
in the 1980s, and decreasing in the 1990s. 


Table 0-2. Flux Boundary Conditions for W (in 109 kg/yr) 



Flux 

Year 

Flux 

Year(s) 

Flux 

l 

5.0 

11 

11.79 

21 

10.61 

2 

5.5 

12 

11.79 

22 

9.43 

3 

6.05 

13 

11.79 

23 

8.25 

4 

6.66 

14 

11.79 

24 

7.07 

5 

7.32 

15 

11.79 

25 

5.89 

6 

8.05 

16 

11.79 

26 

4.72 

7 

8.86 

17 

11.79 

27 

3.54 

8 

9.74 

18 

11.79 

28 

2.36 

9 

10.72 

19 

11.79 

29 

1.18 

10 

11.79 

20 

11.79 

30-40 

0.0 


We graphed W versus XI in Figures 0-16, 0-17, and 0-18 as a function of year for AER, 
GSFC, and WASH, respectively. We wanted to find out how many years were necessary in a 
model simulation before the distributions of a tracer input as a flux (W) compared well with the 
distributions of a tracer at steady state (XI). Since the flux of W varied temporally, its values had 
to be normalized to maximum values of 1 .0 ppbv in order to ease the comparison. 




The points of W versus XI did not line up on the straight line until year 30 in all three models. 
There is a wider extent of points from the WASH model for year 30 (Figure 0-1 8f). This is a 
result of the use of a repeating 4-year cycle of transport in the WASH W simulation (NMC 
temperatures taken from 1979, 1980, 1981, and 1982) versus a fixed 1 year used in the WASH 
XI simulation (NMC temperatures taken from 1980 only). 


These results of W versus X 1 indicate that if models represent the transport in the atmosphere 
in a reasonable way, then chlorofluorocarbons might be expected to take up to 30 years or longer 
from the time they are first released to correlate well with other long-lived constituents (such as 
N2O). 


CONCLUSIONS 

Idealized tracers XI and X2 help in isolating the effect of transport on lifetimes and tracer 
distributions. Models predicting short lifetimes for XI and X2 (NOCAR, GSFC, CALJPL) 
transport constituents more vigorously upward, while models predicting long lifetimes for XI and 
X2 (WASH, DUPONT, AER) transport constituents much more sluggishly. The range of 
lifetimes for XI is from 92 to 127 years and of X2 is from 118 to 161 years. These ranges 
illustrate the approximate bounds of present-day two-dimensional models for simulating these 
idealized tracers. The ratio of the largest to the smallest lifetime for XI is 1.38 and for X2 is 1.36. 
Generally, a short lifetime for XI implies a short lifetime for X2 when compared with other 
models. Each model transports constituents in its own individual way: The latitude-z* 
distributions and temporal behavior of XI and X2 are different from each model. These model 
simulations have indicated the large variability in model predictions from only their transport 
differences. 
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FIGURE CAPTIONS 


Figures 0-1 to 0-10. X2 (in ppbv, represented by large dashed lines) in a) March, b) June, c) 
September and d) December for (0-1) AER, (0-2) CALJPL, (0-3) CAMED, (0-4) DUPONT, 
(0-5) GSFC, (0-6) ITALY, (0-7) LLNL, (0-8) NCAR, (6-9) NOCAR, and (0-10) WASH 
models. Contour levels are l.E-5, l.E-4, 0.001, 0.01, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, and 
0 9 ppbv. Vertical velocities (w in mm/sec) are also plotted; solid lines indicate positive velocities 
and dotted lines indicate negative velocities. For contour values in the w field see Figure A-8. 


Figure 0-11. Tracer XI (ppbv) from all 10 models as a function of month for (a) equator and 
20 km, (b) equator and 30 km, (c) 35N and 20 km, (d) 35N and 30 km, (e) 70N and 20 km, and 
(f) 70N and 30 km. 

Figure 0-12. Tracer X2 (ppbv) from all 10 models as a function of month for (a) equator and 
20 km (b) equator and 30 km, (c) equator and 40 km, (d) 35N and 20 km, (e) 35N and 30 km, (f) 
35N and 40 km, (g) 70N and 20 km, (h) 70N and 30 km, and (i) 70N and 40 km. 


Figure 0-13. Tracer XI (ppbv) from all 10 models in March as a function of altitude for (a) 
70S, (b) 35S, (c) the equator, (d) 35N, and (e) 70N. 

Figure 0-14. Tracer X2 (ppbv) from all 10 models in March as a function of altitude for (a) 
70S, (b) 35S, (c) the equator, (d) 35N, and (e) 70N. 


Figure 0-15. Tracer XI versus tracer X2 for March from (a) AER, (b) CALJPL, (c) CAMED, 
tdl DUPONT (e) GSFC, (f) ITALY, (g) LLNL, (h) NCAR, (i) NOCAR, and 0) WASH. Model 
results are represented by circles (90S?40S), squares (30S-30N), and triangles (40N-90N). The 
straight line indicates the ratio of the lifetime of XI to X2 for each model. 


Figure 0-16. Tracer W versus tracer XI for March from AER for years (a) 5, (b) 10, (c) 15, (d) 
20 (e) 25 (f) 30 (g) 35, and (h) 40. Model results are represented by circles (90S-40S), squares 
(30S-30N), and triangles (40N-90N). The straight line indicates the ratio of 1.0 since W and XI 
have the same loss frequency. 


Figure 0-17. Tracer W versus tracer XI for March from GSFC for years (a) 5, (b) 10, (c) 15, 
(d) 20 (e) 25, (0 30, (g) 35, and (h) 40. Model results are represented by circles (90S-40S) 
squares (30S-30N), and triangles (40N-90N). The straight line indicates the ratio of 1.0 since W 
and XI have the same loss frequency. 


Figure 0-18. Tracer W versus tracer XI for March from WASH for years (a) 5,(b) 

(d) 20 (e) 25 (f) 30 (g) 35, and (h) 40. Model results are represented by circles (90S-40S), 
squares (30S-30N), and triangles (40N-90N). The straight line indicates the ratio of 1.0 since W 
and XI have the same loss frequency. 
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Figure 0-1 
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Figure 0-15 (cont.) 
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